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Abstract 

We study CP violation and the contribution of the strong kaon-pion interactions in the three body 
B — > i^TT+TT" decays. We extend our recent work on the effect of the two-pion S- and P-wave interactions 
to that of the corresponding kaon-pion ones. The weak amphtudes have a first term derived in QCD 
factorization and a second one as a phenomenological contribution added to the QCD penguin amphtudes. 
The effective QCD coefficients include the leading order contributions plus next-to-leading order vertex and 
penguins corrections. The matrix elements of the transition to the vacuum of the kaon-pion pairs, appearing 
naturally in the factorization formulation, are described by the strange Ktt scalar (S'-wave) and vector 
(P-wave) form factors. These are determined from Muskhelishvili-Omnes coupled channel equations using 
experimental kaon-pion T-matrix elements, together with chiral symmetry and asymptotic QCD constraints. 
Prom the scalar form factor study, the modulus of the iV'o(1430) decay constant is found to be (32 ±5) MeV. 
The additional phenomenological amplitudes are fitted to reproduce the K-k effective mass and helicity 
angle distributions, the B — > i^*(892)7r branching ratios and the CP asymmetries of the recent data from 
Belle and BaBar collaborations. We use also the new measurement by the BaBar group of the phase 
difference between the P° and 5° decay amplitudes to ii'*(892)7r. Our predicted B"^ Kq{1A'^Q)'k'^ , 
ii'Q(1430) K'^-K^ branching fraction, equal to (11.6 ± 0.6) x 10~®, is smaller than the result of the 
analyzes of both collaborations. For the neutral B^ decays, the predicted value is (11.1 ± 0.5) x 10~^. In 
order to reduce the large systematic uncertainties in the experimental determination of the B — > i^'g (1430)7r 
branching fractions, a new parametrization is proposed. It is based on the K-k scalar form factor, well 
constrained by theory and experiments other than those of B decays. 

PACS numbers: 13.25.Hw, 13.75Lb 
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I. INTRODUCTION 



Rare two-body and quasi two-body charmless hadronic decays of B mesons [1] are a rich field 
for tests of the standard model and QCD [2]. Furthermore three-body charmless hadronic B 
decays provide an interesting ground, not only for searches on CP violations but also to study 
hadronic physics [3]. Strong interaction effects, in particular through the presence of two-body 
resonances and their interferences, can influence weak decay observables. Strong interaction phases 
arc necessary for the occurrence of CP violation and it is essential to have a description, as reliable 
as possible, of the interactions between the detected hadrons. Dalitz-plot analyzes allow to extract 
effective mass and angular distributions of the produced meson pairs. If one pair is created in 
two (or more) different states, one can see specific interference effects leading to additional and 
interesting CP asymmetries. These arise from the variations of the strong phases with energy, such 
variations are absent in the two-body B decays where the energy is fixed. The meson-meson final 
state interactions must be addressed using theoretical constraints, such as unitarity, analyticity and 
chiral symmetry, and experimental data from processes other than B decays. Then, for a given 
P-decay, the comparison between the theoretical model and experimental results will determine 
the strong phases needed to generate the measured direct CP asymmetries. 

Recently, BaBar and Belle Collaborations have performed detailed Dalitz plot analyzes for 
different B Kit^tt~ decays [4-12]. One observes an accumulation of events for tttt or K-k 
effective masses lower than 2 GeV and in particular the presence of the scalar mesons /o(980), 
if^(1430) (and to a lesser extent /o(1500)) and of the vector mesons /9(770)°, if*(892). The event 
distributions of the Dalitz plots are usually studied using the isobar model in which the decay 
amplitudes are parametrized by sums of Breit-Wigner terms and a background. In Refs. [4, 7, 10- 
12], an effective range nonresonant component has been used in the Kti 5- wave amplitude. The 
aim was to ameliorate the description of the low Kir effective mass spectrum. 

An important breakthrough in the theory of B decays recently achieved is the confirmation of 
the validity of factorization, as a leading order approximation in an expansion in inverse powers of 
the h quark mass mi, [13]. This concerns the B decays into two mesons and was later reformulated 
using the soft coUinear effective theory approach to QCD [14]. Detailed comparisons between 
theory, based on QCD factorization (QCDF), and experiment were made in the case of decays 
into two pseudo-scalar mesons and one pseudo-scalar and one vector meson [13, 15]. Agreement is 
generally quite fair. However, some phenomenological parameters are introduced. 

In this paper, we study the decays of the B into three mesons, B —>■ i^Tr+Tr", for which, to 
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our knowledge, no proof of factorization has been given. However, we restrict ourselves to specific 
kinematical configurations in which the three mesons are quasi aligned in the rest frame of the 
B. This condition is met, in particular, in the low effective Kir mass region (< 2 GeV) of the 
Dalitz plot where most of the resonant structures are seen. We will denote such processes as 
B {Ktt)tt where the mesons of the i^vr pair move, more or less, in the same direction. Three- 
body interactions arc expected to be suppressed in such conditions. Then, it seems reasonable to 
assume the validity of factorization for this quasi two-body B decay [16] where we assume that the 
(Ktt) pair originates from a quark-antiquark state. 

In a previous work [17], the decays B — (7r"''7r~)5 K as well as B ^ {KK)s K, where the two 
comoving mesons in the (7r'*"7r~) and (KK) pairs are in S-states, were studied using an approximate 
construction of relevant scalar form factors proposed in Ref. [18]. The decays B (7r+7r~)p K, 
with the two pions in a P state, were subsequently studied in Ref. [19]. Here, we focus on the 
decays B^ {K^it^)tt^ , B^ {K^^it^)tt^ and B^ {K^tt^)'k^ . In the factorization approach 
the amplitudes can be expressed as the product of effective QCD coefficients [13] and the two matrix 
elements of the vector currents {KTr\q^'^ s\0) and {Tr\q^i,b\B) , with q = u ot d. One conspicuous 
consequence is that the Ktt pair is restricted to be in either an S or a P state; no D or higher 
waves being allowed. This fact is indeed supported by the Belle experiment of Ref. [9]. In their 
analysis on P° — K^tt+tt" decays, one clearly observes a vector if* (892)+ and a scalar ifQ(1430) 
resonances, the signal for the tensor if|(1430) being small. 

The matrix element of the vector current {KTT\q^'^s\0) involves two functions of the Ktt effective 
mass squared, the strange scalar and vector form factors. We will perform a construction based 
on general properties of analyticity, QCD asymptotic counting rules [20] and using accurate ex- 
perimental data on Ktt scattering, both in the elastic and inelastic region. Such data have been 
obtained in the relatively high statistics LASS experiment on K~ proton interaction [21-23]. A 
specific feature of Ktt scattering at medium energies is that inelasticity is dominated by two-body 
or quasi two-body channels. More specifically, it was shown by LASS that inelasticity in the S- 
wave is dominated by the Krj' state and in the P-wave by the if*7r and Kp states. Combining 
dispersion relations with unitarity relations then leads to a set of coupled integral equations for 
the form factors. These equations constitute a generalization of Watson's theorem of final state 
interactions, well known for a single channel case [24]. 

An analogous system of equations was studied and solved for the first time by Donoghue, Gasser 
and Leutwyler [25] in the case of the pion scalar form factor. More recently, this method was 
applied to the Ktt scalar form factor in Ref. [26]. For our purpose, we have redone the calculation 



4 



of Ref. [26] and extended this framework to the case of the vector form factor. Both of these form 
factors, needed in our work on B ^ {Ktt) tt decays, are treated exactly in the same way. This 
provides an unified treatment of different B decays through coupled final-state channels and also 
of r ^ KttUt decays [27, 28] using form factors constrained by accurate results on Ktt scattering. 
It could allow, for instance, to give predictions for B [Kr]')s vr, {K*tt)p vr and {Kp)p vr decays 
as was done in Ref. [17] for B {KK)s K decays connected to i? — > (7r"'"7r~)5' K decays. 

The introduction of form factors, constrained by theory and other experiments than B decays, 
is an alternative to the use of the isobar model. This latter approximation violates unitarity and 
the information about resonances, present in the final states, can be distorted by other nearby 
resonances due to interference. In our approach, we use the complex pole definition of a resonance 
which allows us to obtain its branching ratio and its decay constant. An alternative way, examined 
here too, is to integrate over the effective mass range where the resonance dominates. Recently 
charmless three-body decays of B mesons have been extensively studied by Cheng, Chua and Soni 
in the factorization scheme [29]. Their calculation proceeds via quasi two-body decays involving 
resonant states and nonresonant contributions. Breit-Wigner expressions are used to describe the 
appropriate resonance effects in the scalar and vector matrix elements. Their K*{892)-jt, p(770)°if 
and KQ(1430)7r branching ratios are too small compared to the data by a factor varying between 
2 and 5. This is also the case in our QCD factorization approach. To improve agreement with 
experiment we introduce phenomenological corrections to the QCD penguin amplitudes. The latter 
could represent (in part) the contribution (not studied here) of the weak annihilation and hard- 
spectator contributions together with their phenomenological components [13]. They could also 
partially come from long distance charming penguin amplitudes [14, 30] which, themselves, could 
arise from intermediate states, reminding that B — > D^*^ D^*^ branching fractions are 

large. Finally, they can come from unknown effects of new physics which could appear in the 6 — ^ s 
quark loop transitions (see for instance Ref. [31]). 

Our paper is organized as follows. In Sec. II we derive, in the QCDF framework, the B 
{K-ir)s,p 7r='= decay amplitudes for 5=*=, B° and 5°. The charged B {b ^ sdd transition for B~) 
decay amplitudes have only penguin diagram contributions, while the neutral {h suu transition 
for B^) have an additional tree diagram. In the weak amplitudes, we include penguin-correction 
terms represented by four complex parameters. The amplitudes are expressed in terms of the 
product of the effective QCD coefficients by the B to tt transition form factor and the K-k strange 
form factors. Section III tabulates the values of the process and scale dependent effective QCD 
coefficients [13] we use in our amplitudes. To the leading-order contribution in a^, we add the 



5 



next-to-leading order short-distance vertex and penguin corrections. Their calculation is outlined 
in Appendix A where we give also their values. 

In Sec. IV we specify the model we use for the B to tt transition form factors. The unitary 
equations satisfied by the scalar and vector strange form factors are also presented. We discuss 
briefly the low energy constraints and the two-channel description {Kn, Kr]') for the S-wave Kn 
scattering. In the case of the P-wavc, the necessary three channel (Kir, K*tt, Kp) description 
is then described. We give the results for the strange scalar and vector Kit form factors. Using 
the complex pole definition of a resonance, a simple and unambiguous separation between the 
background and the resonant contributions of the S-wave and P-wave amplitudes, is given. This 
allows us to determine in a unambiguous way the B iCQ(1430)7r and B K*(892))7r branching 
fractions. A detailed derivation of the S-wave T-matrix elements, necessary to calculate the strange 
scalar form factor, is presented in the Appendix B. We calculate the values of the decay constants of 
i^*(1430) and i^Q (892) from the knowledge of the vector and scalar Ktt form factors in Appendix C. 
In Appendix D, we show how the two body amplitudes B i^Q(1430)-7r and B i^*(892))7r are 
related to the three-body ones. Effective decay constants for iir*(1430) and Kq{892) are also 
calculated. 

In Sec. V we describe our fitting procedure on the additional penguin parameters. We further- 
more compare the results of our fit to the experimental mKn mass and helicity-angle distributions. 
We present our fitted values for the branching ratios and CP asymmetries for the B — ^ K*(892)7r 
and B i^Q(1430)-7r decays. Note that, in the case of the scalar meson production, our branching 
ratios are predictions. Discussion of the results and comparison with experimental analyzes are 
also given. A summary and some conclusions are presented in Sec. VI. 

II. DECAY AMPLITUDES AND PHYSICAL OBSERVABLES 

The amplitudes for the non-leptonic decays of the B meson are given as matrix elements of the 
effective weak Hamiltonian 

^^ff = -^Y. + <^20f + 5Z CiOi + C7^07^ + CsgOsg\ + h.c, (1) 

^ p=u,c 1=3 

where in the case of strangeness S = ±1 final states 

A„ = V^bV:,, Xc = V,bV:„ (2) 

the Vppi being Cabibbo-Kobayashi-Maskawa quark-mixing matrix elements. For the Fermi coupling 
constant Gp we take the value 1.16637 10~^ GeV"^. In this work we use A„ = 3.55 x 10~^-i 7.49 x 
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10~^ and Ac = 4.05 x 10~^ + i 6.5 x 10"'^. The Ci(/x) are the Wilson coefficients of the respective 
four-quark operators Oi{fx) at a given renormahzation scale /i. The explicit expression of the 
operators Oi may be found e.g. in Ref. [32] . Studies of B decays into two-body [32] and quasi- two- 
body [13, 33] final states have been performed in the QCDF framework. These studies show that 

naive factorization is a useful first order approximation which receives corrections proportional to 
the strong coupling constant 0^(7715), as{y^ AgcD'mb) and in inverse powers of nib [2]. Here, we 
perform a heuristic extension of these results to a class of three body decays B — > {K'!r)7r. 



A. Charged B decay amplitudes 

We focus on the process B~ — {K~Tr'^) tt". To illustrate our approach, let us write the matrix 
elements of the penguin operators O3 and O4 at leading order factorization^ 

(7r-(K-7r+)|C303 + C4O4I5-) = a4{7r-\dj''{l - j^)b\B-){K-Tr+\sj,{l - 75)c?|0) (3) 

with 

a4 = C4(/x) + i- C3(/i), (4) 

where iVc = 3 is the number of colors. In this approximation, the dependence of the amplitude as 
a function of the two Dalitz-plot variables rriKTn "^tttt is completely determined in terms of K — tt 
and B — ^ TT form factors. We will probe this prediction by employing a careful determination of 
the Ktt form factors described in detail in Sect. IV and in Appendix B. Our main assumption 
will be that the corrections to naive factorization can be absorbed into effective-mass independent 
modifications of the parameters a,. We will borrow parts of these corrections from quasi-two-body 
calculations and also append a phcnomenological part. 
The B to TT transition matrix element is written as 



{7r-{p,-)\dY{l-l')b\B-{pB-)) 



Ml -ml "-^^ _2 



fr--iq^) + MBj!^g.fB-.-^g2^^ (5) 



where f^ i {q^) are the scalar and vector B to tt form factors. The four-momentum transfer 
is 

q = PB- -P-K- =Pk- (6) 



^ In this derivation we have assumed that the {K tt"*") pair originates from a quark-antiquark state. 
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with Ptt-) Ptv+i Pk- ^'iid pb- being the four momenta of the negative and positive pions, of the 
K~ and of the B~ mesons, respectively. In an analogous way, the matrix element for the transition 
from vacuum to the K~Tr'^ state is given in terms of the scalar /(f- '^^ (q^) and vector ff^ '^^ {q^) 
form factors by 



{K-{pK-)7T+{p,+ )\s-f,{l - 75)d|0) 



{PK- -P'K+)v - 



Above, Mb, rrix and rriT^ are the masses of the charged B mesons, kaons and pions, respectively. 
From Eqs. (5) and (7) one obtains 

{7r-\dY (1 - 7') ^'l^") {K-7T+\sj,{l - 75)d|0) 



21 



m 



K-n- 



{Ml -ml] 



, (8) 



where mK--n- and m^+T^- are the K n and tt+tt effective masses. Note that q"^ is the square of 
the K~Tr~^ effective mass. In the K~Tr~^ center of mass system 



m 



K-TT- 



{Ml - ml) 



mj^-m- 



= 4p^+ • Ptt- = 4|P7r+ 1 IPtt- I cos 9, 



(9) 



where 9 is the angle between the tt'*' and 7r~ three-momenta. The first term of the right hand 
side of Eq. (8) corresponds to the K'tt'^ S'-wave contribution while the second term to that of the 
P-wave since it depends linearly on cos 9. One observes that there are no contributions from I > 2 
partial waves. 

Finally, wc introduce the matrix element of the complete weak effective Lagrangian. We write 
the Kir S-waye contribution in the following form 



Gf. 



Ms ^ {n- iK-7r+)s\H,ff\B-) = ^{Ml-m^^ 



X ( al{S) - + ci] + Ae ( al{S) - + 0% 



2g2 



(mfe - md){ms - ma) 



(10) 



Here the kinematical q^ dependence associated with the terms 04 — a^g/^ + c^ {p = u, c) is different 
from that in front of the Og — + ones. In the latter case, due to a Fierz transformation, 
the matrix elements of scalar rather than vector current are involved. This dependence has an 
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important consequence for the behavior of the effective mass distributions (see Sec. VB 1). For the 
P-wave one has, 



Mp p,- • p,+ = (tt- {K-7r+)p\Heff\B-) = 2V2Gf p,- ■ p,+ /f ^(g') 



K (<(P) - ^ + c^) + A. (<(P) - ^ + c^) 



(11) 



rub fv 

In Eqs. (10) and (11) a^{S) and a^(P) are the leading order factorization coefficients to which 0{as) 
vertex and penguin corrections are added in the quasi-two-body approximation of pseudoscalar- 
scalar, PS, or pseudoscalar- vector, PV, final states (see Sect. III). These coefficients will be 
evaluated at the scale n = m,},. The term proportional to fyilA/fv has been inferred from a 
similar term which arises at order ag in the B PV amplitudes [13]. In the following, we take 
the i^*(892) decay constant fv = 218 MeV and /^(rrife) = 175 MeV with mK* = 893.8 MeV. 
We have further assumed that the corresponding P-wave form factor (which is associated with the 
antisymmetric tensor current) is simply proportional to the one associated with the vector current. 

It was observed in Ref. [13] that the calculated 0{as) and l/m;, corrections are insufficient to 
explain the experimental branching rate for the B K*Tr decay. We have therefore allowed for 
four additional complex terms c\, 0%, Cg and Cg which one could partly interpret as non-perturbative 
contributions in the penguin diagrams. The other part of these coefficients could represent hard 
spectator interaction in the pertubativc regime and annihilation terms. Annihilation diagrams, due 
to divergences inherent to the form of the twist-three amplitudes employed in QCDF [13, 32], are 
commonly parametrized by complex amplitudes very similar to these QCD penguin corrections. 
This complex parametrization is also the case for the end-point divergence of the hard-spectator 
amplitudes. The above parameters will be determined by performing detailed fits to the Dalitz 
plot in the region m^Tr < 1-8 GeV, mTTTr > 1-5 GeV where one believes that the QCDF formalism 
could apply. 

The complete amplitude for the B~ {K~Tr'^)Tr'~ decay is 

M- =Ms+MpP^--p^+. (12) 
Charge conjugation of M.~ gives the B~^ —> (i^+7r~)7r+ decay amplitude 

7W+ = M-{Xu a;, Ac ^ A*, 5" ^ B+, K- K+, tt"^). (13) 

Let us specify here the value of the meson masses used in this work: Mb = 5279.2 MeV, 
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rriK = 495.66 MeV (averaged and masses) and = 139.57 MeV. For the quark masses, 
evaluated at the scale fi = mj,, we take = 4.2 GeV, = 84 MeV and m„ = = 3.4 MeV. 

B. Neutral B decay amplitudes 

In the 5° — > {K^^7T^)s,p 7r+ decay we have the quark transition b — suu. The derivation is 
similar to that just described above for the B" (i^~7r+)s,p vr~ case. However, there is a tree 
diagram ai contribution. One obtains for the ^-wave amplitude 

2 2 

■^S = \^ ^ )s\-tieff\-D } - —i=[Mgo - mj 2 JO W ) JQ KQ ) 

y 2 q 



X { A„(ai + aliS) + a^S) + 4) + Xcial{S) + aloiS) + 4) 
2g2 



(mft - mu){ms - my) 
For the P-wave amplitude one has 



[A„K(5) + al{S) + cl) + X,{al{S) + a^(5) + c^)] . (14) 



X |a„ (ai + <(P) + <o(P) + c|) + Ae {al{P) + a5o(P) + 4) 

+ 2!!!^^:^ [A„ «(P) + al{P) + ce") + A^ (a^(P) + ^(P) + c^)] |. (15) 
The full amplitude for the — > (^'^7r~)7r+ decays is 

>(' = A4^ + >(PP.--P.+ - (16) 
The charge conjugation of the MP amplitude gives the P° — > (iir°7r"'")7r~ decay amplitude as 

= M^{Xu ^ a;, Ac ^ A*, po ^ po, ^ i^o, TT^ ^ 7r±). (17) 
Isospin symmetry leads to 

/(fi '^^ {l') = /cfi (^^) = fo^,r~ il') = foT" (18) 

and to 

fo^'i^) = foT"(^') = /oT^'(«') = foT'i^')- (19) 

In Eqs. (14) and (15) we have introduced the same phenomenological parameters as those for the 
charged P-decay amplitudes, Eqs. (10) and (11). 
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C. Physical observables 



The density distribution of the Dalitz plot for the B" K~7r+7r~ decay^ can be expressed in 
terms of the K~7r+ and tt+tt" effective masses, the latter being related to cos^ as seen in Eq. (9). 
The double differential B~ K~Tr'^Tr~ decay rate reads 



mK-TT+lPn+l IPtt- 



■\M- 



dcos9dmK-n+ ~ 8(27r)3M| 
where M~is given by Eq. (12). The moduli of the vr"'" and tt" momenta are 

1 



2m 



m 



K-TT+ 



{triK + m^) 



m 



(20) 



(21) 



1 



IPtt- 



(22) 



2m^-^4 

In the experimental analyzes, the helicity angle 9h is usually defined as Oh = — 0- Integrating 
the double differential distribution of Eq. (20) over cos 6 gives for the differential effective mass 
branching fraction 



dB- 



+-^\Pn+\ IPtt-I \Mp\ 



dmK-^+ Tb- 4(27r)3M3 
where Tg- denotes the total width of This is a sum of the S- and P-wave contributions 



(23) 



dB- 



dB: 



+ 



dB\ 



dmK-n+ dmK-n+ dmK-n+ 
The CP violating asymmetry for the charged B decays is defined as 

B- - B+ 

Acp = 



(24) 



B-+B+ 

The integration of the double differential decay rate over mK--K+ within the range (mmim rn 
gives the angular distribution 

dB- 



(25) 



dcosO 



A + BcosO + C cos^ 



(26) 



where 



A 



J m 



dm 



mK-n+lPTT+l IPtt-I I . . 



8(27r)3M| 



2 

51 ' 



(27) 



^ Note that here and in the following, for simplicity, we suppress the parentheses between the two first mesons of 
the three-meson final state. 
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dmK--jv+ 



mK-n+\Pn+f \Pn- 

8(27r)3M| 



\Mp 



(29) 



Similarly one can derive the above observables given in Eqs. (20), (23) and (26) for the 5+ 
K+TT-'!r+, S° ^°7r-7r+ and S° K^tt+tt- reactions. 

D. B ^ K*(1430)7r and B K*(892))7r amplitudes 

The scalar Kq{1A30) and the vector K*(892) resonances are quite visible in the experimental 
effective Ktt mass distributions of the B — > Ktttt decays [4—12]. Branching fractions for the 
B —> KQ(1430)7r and B —> K*(892)7r reactions have been extracted by the experimental groups 
within the framework of the isobar model. Here, we will make use of the complex pole definition of 
a resonance in scattering theory. This method allows one to perform in a simple way a separation 
between the resonance contribution (defined to correspond to the pole part) and the background 
contribution in the amplitude. 

We obtain the B K^{U30)tt amplitudes by replacing in Eqs. (10) and (14) /^""(q^) by 
/o°'^(5^) defined in Eqs. (45), (46) and (47). Similarly the B K*{892)tt amplitudes are given by 
replacing in Eqs. (11) and (15) ff^^iq"^) by fl°^'^{q^) defined in Eqs. (62) and (63). Integration of 
the differential branching fraction (23) over the effective mass mxTi of the pole parts of the scalar 
or vector form factors will determine the B~ KQ'^(1430)7r~ and B~ K*'^(892)7r~ branching 
fractions, respectively. In this approach, we can also determine the decay constants of the vector 
and the scalar resonances: this is exposed in Appendix C. 

III. EFFECTIVE QCD AMPLITUDES 

The short-distance physics of the weak decay amplitudes b sdd or 6 — > suu is codified in 
the effective QCD amplitudes a"'^(;[i), each of which corresponds to a particular decay topology. 
The available QCDF calculations (see for instance Refs. [13, 32, 33]) concern two-body final states. 
They apply to our case of a three-body final state (K7r)7r when the effective mass of the K-k sub- 
system coincides with the mass of a vector or a scalar resonance. We will make the approximation 
to use these calculated corrections also away from the resonance masses. In this work, we take into 
account only vertex and penguin corrections. At the leading order, these amplitudes are universal 
(i.e. do not depend on quark flavor) and given by 
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TABLE I: Leading and next-to-leading order effective QCD amplitudes ai{mb) and a^(m(,) for the JsTq (1430)7r 
and /i:*(892)7r final states [see Eqs. (30) and (31)]. 





S - 


■^i('o*( 1430)77 


B - 


-> i4:*(892)7r 




a^TOb) 


of (rrib) 


ai{mb) 




fll 


1.018 


1.029 + i 0.063 


1.018 


1.045 + i 0.014 


al 


-0.031 


-0.061 0.023 


-0.031 


-0.030 0.015 




-0.031 


-0.069 + i 0.057 


-0.031 


-0.035 - i 0.005 


"6 


-0.039 


-0.042 0.014 





-0.006 - i 0.002 


ag 


-0.039 


-0.045 - i 0.004 





0.002 + i 0.009 


"8 


0.00044 


0.0005 -i 0.0001 





-0.0001 + i 0.0 


ag 


0.00044 


0.0005 0.0 





-0.0 + i 0.0001 


ti 

"lO 


-0.0015 


-0.002 + i 0.003 


-0.0015 


0.0001 + i 0.0006 


«10 


-0.0015 


-0.002 + i 0.004 


-0.0015 


0.0001 + i 0.0007 



aiifJ') 



Ni{M), (30) 



where the plus (minus) sign corresponds to odd (even) values of the index i. One finds Ni{M) = 
if the emitted meson M (either the K*{S'd2) or the ii'g(1430), which do not include the spectator 
quark) is a vector one and if i = 6,8, otherwise Ni{M) = 1. The next-to-leading order vertex and 
penguin corrections are calculated following Ref. [33] for i^Q(1430)7r final states and Ref. [13] for 
if*(892)7r ones. We write 

<(/.) = a,(/x) + + PiiM), (31) 

with Cf = (A^c ~ l)/2A'c. The relevant Wilson Cj(;u) coefficients were calculated in next-to-leading 
order logarithmic approximation and we take their values from Table 1 of Ref. [32] at the natural 
characteristic scale fj, = mi,. The second and third terms of Eq. (31) refer to the vertex and penguin 
corrections and for completeness we describe in Appendix A the main steps of their calculations. 
In Table I we list the values of the coefficients a^{mb) with i = 1,4, 6, 8, 10, used in the weak decay 
amplitudes discussed in the present paper. These have been obtained adding to the aj(mb) values 
(given in this table) those of the vertex and penguin corrections listed in Table VI. The next-to- 
leading order corrections are relatively small except in the case of the B irQ(1430)7r decays for 
04 and a|. These arise mainly from the penguin terms. 
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IV. SCALAR AND VECTOR Btt AND Ktt FORM FACTORS 



A. Bit form factors 

The scalar and vector Bit transition form factors fo^iq^) and fi^{q^) were introduced in Eq. (5) 
and are one of the ingredients in QCDF. These form factors contain parts of the non-perturbative 
physics which stems from the hadronization of quark currents. In the appUcations considered here, 
the variable is timelike and remains small compared to the physical threshold (me + m.T^)'^. 
The form factors in this region are real and we expect their variation to be slow. Several studies 
dedicated to these transition form factors have been carried out; on the phenomenological side, we 
refer to the various approaches to heavy-to- light transition amplitudes. These include light-cone 
sum rules [34], light-front [35], simple non-relativistic [36] or relativistic quark models [37-40]. 
Ab initio approaches such as lattice-regularized QCD have also been used to determine Btt form 
factors, however the lattice results are for large momenta > 10 GeV^ and need to be extrapolated 
to low momenta [41, 42] by means of pole dominance models. 

More recently, a comprehensive set of B-meson heavy- to light-transition form factors, cal- 
culated with truncated transition amplitudes based on Dyson- Schwinger equations in QCD, was 
reported in Ref. [43]. The methods of Refs. [34, 35] yield form factors for a small domain of 
time-like momentum transfer q'^, while those in Refs. [37-39] apply to the entire range of physical 
momenta. However, this is only possible by employing functional extrapolations of the transition 
form factors foJ{q^ = 0) to the time-like g^-range. On the other hand, the form factors obtained 
from double dispersion relations of spectral densities in the relativistic quark model [40] or us- 
ing Dyson-Schwinger equations in QCD [43] are calculated for all physical q^ values. A typical 
value found in [40] is fi'"{q^ = 0) 2± 0.2 which is in agreement with the Dyson-Schwinger result 
fi^iq^ = 0) = 0.24 as well as with lattice data extrapolations [42]. 

In the above studies, the dependence of the Btt transition form factors, as varies from 
threshold up to about (1.8 GeV)^, is found to be small. This, of course, is in contrast with the case 
of the Kir scalar or vector form factors which we discuss below. In practice, we take the following 
constant values: f^^q^) = /(fn^K5(i430)) = 0-266 and /f'^(g2) = /f '^(m^.^^gg^)) = 0-250. 

B. Ktv form factors 

Another important ingredient of the QCD-factorized B-decay amplitudes are the Ktt scalar 
and vector form factors. These also appear in semileptonic decays like r Kirvj- or K ^ ttIui. 



14 



Use of analyticity and unitarity allows one to relate them to the S and the P wave Kir scattering 
amplitudes. Indeed, rather accurate information on Kw scattering is available, in particular from 
the high statistics experiment of Estabrooks et al. [44] and from the LASS Collaboration [23] . This 
approach served to determine the scalar form factors of the pion and kaon in Ref. [25]. Jamin, 
Oiler and Pich applied it recently to evaluate the Kit scalar form factor [26]. We perform also its 
calculation following, essentially, their work. We then develop an extension of this construction to 
the vector case, so both form factors are handled in a similar way. The main approximation is to 
reduce the sum over the inelastic many-body channels in the unitarity equations to a finite number 
of two-body channels. This is supported by the experiments performed by the LASS Collaboration 
in the center of mass energy range m^Tr < 2.5 GeV [21-23, 45]. They find that the inelasticity in 
the S'-wave is dominated by one channel, Krj', and in the P-wave by two channels, K*Tr and Kp. 
Then, from fits to the experimental data, one constructs a 2 x 2 scattering T-matrix for the S wave 
and a 3 X 3 one for the P wave. The form factors satisfy a set of n coupled, homogeneous singular 
integral equations with a kernel linear in the T-matrix (n = 2 for the S wave and n = 3 for the P 
wave). The mathematical properties of such equations are derived in Muskhelishvili's book [46]. 
In particular the number of independent solutions N is given by the index of the integral operator 
which can be expressed in terms of the sum of the S-matrix eigenphases dj{t) of the S'-matrix, 

n 

^[5,(oo)-<5,(0)]=Ar7r. (32) 

N is also the number of independent conditions that one must impose on the form factors in order 
to determine them from the integral equations. We will use conditions at t = or near t = 
derived from chiral symmetry. Asymptotic conditions will also be used. 

1. Scalar form factor 

Following the experimental analyzes, we assume the dominance of a single inelastic channel 
Krj' ^. The Kri' matrix element of the vector current can be written in terms of two form factors 

/^^'''(t)and/5^'''(t)as 

{K+\urs\v') = ff'^'mpK+Pr,')^ + /^^'''(i) W -P^')^ (33) 

^ The authors of Ref. [26] have studied the influence of an additional Kr] inelastic channel and found that it is rather 
small. 



15 



3.5 
3 

2.5 
2 

1.5 
1 

0.5 





/a-//. = 1-203 

./a-//^ = 1-193 / '. 

IkIL = 1-183 /■A\ 




250 
200 
I 150 
100 
50 




JkIL = 1-203 
JkIL = 1-193 
/a/A = 1-183 



0.5 1 1.5 



2.5 



0.5 




^ft ( GeV) 



yft ( GeV) 



FIG. 1: Modulus (left) and phase (right) of the strange scalar form factor '^^(0 obtained by solving a 
two-channel Muskhelishvili-Omnes equation system. Variation with the input Fi{/S.ktt) = /k/Ztt — 3.1 x 10^"^ 
at the Cheng-Dashen point A^tt = "rn]^ ~ 'ti^ is illustrated. In the present work we use the form factor 
corresponding to /aV/tt = 1.193. 

where t = {pK—Pri')'^ = q'^- The coupled channel integral equations will involve the two components 



F,{t)^V2f^'--\t) = fr-\t), 



2 2 

\f^ ' it) + 



m\ — rn^ 



(34) 



We have used here isospin symmetry to express the form factor '^^(i), introduced in Eq. (7), 
in terms of f^^'^^ (t). 

The form factors Fi{t) and F2{t) are analytic functions in the complex t plane with a cut along 
the real axis (m,r + itt-k)^ 1^ t < oo. Asymptotic counting rules in QCD [20] imply that the 
dispersion relations satisfied by the functions Fi{t) converge without subtractions 



1 

m) = - / 



Im Fi{t')dt' 



i = 1,2 . 



(35) 



Unitary equations of the two coupled channels Ktt (i = 1) and Kr]' (i = 2) allow one to express 
the imaginary parts as follows (see [25, 26]) 



ImFi(t) = e{t - ti)^^^ TA(t) Fi(t) + 9{t - t,)'^-^^TUt) F^it), 



ImF2(t) = eit- t,)^^^ Ti*2(t) Fi(t) + 9{t - t2) 



n2{t)F2{t), 



(36) 



where 



h = {'rriK + "t-tt)^, *2 = {rriK + m^iY 



(37) 
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and qKn{t), qKr)'{t) are the center of mass momenta. The set of integral equations obtained 
by inserting Eqs. (36) into (35) are often called Muskhelishvili-Omnes equations. The detailed 
determination of the matrix elements Tn, Tu, T22 is given in Appendix B. We use experimental 
data and theoretical constraints at low energy, in particular we employ a systematic combination 
of the chiral and 1/Nc expansions [47, 48]. At medium and high energies, we employ a X-matrix 
parametrization as in Rcf. [26] which guarantees unitarity. This provides a smooth interpolation 
with index A'" = 2 in the asymptotic region, where no experimental data exist [see Eq. (32)]. One 
must therefore provide two initial conditions for the form factors. As in Ref. [26], we use the values 
of Fi {t) at t = and at the Cheng-Dashen point t = — which are precisely constrained by 
chiral symmetry 



Fi(0)= 0.961, Fl(m|f-m^) = ^-3.1•10-^ (38) 

where fx and /tt are the kaon and pion decay constants, respectively. The value at t = is derived 
from chiral perturbation theory at order [49, 50] and includes an estimate for the corrections 
made in Ref. [51]. The value at t = m\ — was obtained in Ref. [50] at order p^. In that case, 
the corrections are expected to be very small, of order 10~^. For the ratio fx/ f-n, the latest 
Review of Particle Physics [52] (see the note by J. Rosner and S. Stone, p. 821) gives 

^ = 1.193 ± 0.002 ± 0.006 ± 0.001 (39) 

and we will use the central value in Eq. (38). We solve the set of integral equations for Fx and F2 
using the numerical method of Ref. [53] . The results for the modulus and for the phase of Fx are 
displayed in Fig. 1. Our numerical results are in a fair agreement with those of Refs. [26] and [27]. 
The phase of the form factor satisfies, as it should, Watson's theorem in the energy region below 
the Krj' threshold. Above this point it displays a sharp drop. Correspondingly, the modulus of 
Fi displays a dip. These structures reflect interference effects between the independent solutions 
which are linearly combined. Due to the linearity of the equations, the dependence on the initial 
conditions can be studied by varying Fi{m\^ — m^) and keeping Fi(0) fixed. The figure illustrates 
the variation of the scalar form factor f^'^ as a function of fx/ f-K- 

2. -K'*(1430) pole part of the scalar form factor 

The separation between resonant and background contributions has some arbitrariness, in 
particular if the background contribution is significant, depending on how each contribution is 
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parametrized. In the case oi B ^ K-mr amplitude, we propose a simple and unambiguous way to 
define this separation based on the well known property that a resonance can be associated with a 
pole of the scattering matrix in the complex energy plane on the second Riemann sheet (e.g. [54]). 
This pole also appears in form factors and current correlation functions. We first study the scalar 
Kq{14:30) resonance and then the vector i^*(892) resonance. 

Let us analyze the scalar form factor fQ^{t). We want to define the extrapolation to the second 
Riemann sheet in the t variable. We recall that we have assumed, in constructing the T matrix, 
that scattering was elastic up to the r]'K threshold. We can write the discontinuity of the form 
factor upon crossing the cut as follows, 

f^^{t + ie) - f^^t - ie) = -2aK^{t + ie)T^^{t + ie)f^%t - ze), (40) 

for t real and in the range (m^ + rriK)'^ < t < (m^/ + m^)^- Here, axnit) = 



{m,K — mT^Y)/t and Tf^ is the 5- wave T-matrix element of the Kir Kir 
process. Tf^(t) satisfies a discontinuity equation similar to Eq. (40), 

Tfi(t + ie) - rfi(t - ie) = -2aK^{t + ie)T^^{t + ie)Tf^{t - ie). (41) 



Eqs. (40) and (41) allow us to find the extension of f^'^ on the second Riemann sheet 

/(f"(t) 
l-2aKnmiiit)' 



/o'(0 = . J'^^JL... , (42) 



which, by definition, must satisfy 

fii{t-ie) = f^^t + ie) (43) 
along the cut. Eq. (42) shows that /o^(t) displays a pole whenever the denominator function 

D{t) = l-2aKAt)Ti\{t) (44) 

displays a zero. In a similar way, one can define the extension to the second sheet of the T-matrix, 
(Tfi)^^ , which because of Eq. (41) has exactly the same denominator function D{t). A point such 
that D{tQ) = corresponds to a pole of (Tf^)^^ and thus can be associated with a resonance [54]. 
Prom this, it is simple to isolate the pole part of the form factor 

where a = dD{t)/dt dX t = to- In the numerator, f^'^{to) can be computed using its dispersive 
representation. 
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Vi (GeV) 

FIG. 2: Modulus of the scalar form factor /J^"'' (t) compared with its pole part, as defined in the text. 

In practice, in order to do so, we must be able to define Tf^{t) for complex values of t. In our 
work we have defined Tf^{t) on the real axis in the range 1.25 <Vi<2.5 GeV from a 2 X 2 K 
matrix fit to the experimental data. By construction, the elements of the i^-matrix have no branch 
cut on the positive real axis. The meromorphic function parametrization which we used should 
be valid in some complex domain of the t variable. It seems reasonable to assume that it remains 
reliable in the region of the -fCg(1430) resonance pole since this pole lies rather close to the real 
axis. Numerically, we find the following result for the location of the pole 

to = (1.9487 - i 0.3825) GeV^ Vt^ = (1.4026 - i 0.1364) GeV. (46) 

These results compare reasonably well with the values of the mass Mr = (1.414 it 0.006) GeV and 
of the half-width Tr/2 = (0.145 ± 0.011) GeV of the K*{U30) given by the PDG [52]. For the 
other quantities needed for /o°''^(t) we obtain 

f^-^{to) = -0.3242 - i 1.4679, a = (0.8381 + i 1.1713) GeV"^. (47) 

We remind that, in the dispersive construction of fQ^{t) from Ref. [26], followed here, the result 
depends on the value of fx/ f-w which controls one of the initial conditions [see Eq. (38)]. The 
number given above for fo{tQ) corresponds to fx/fn = 1-193 which we used as central value. 
Varying fx/ fw we would obtain: 

f^^'ito) = -0.4901 - i 1.6652 for ^ = 1.203, 

Jn 

f^^{tii) = -0.1586 - i 1.2710 for ^ = 1.183. (48) 
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Figure 2 shows the modulus of /q" ^(i) compared with the modulus of {t) and that of the 

"background" which we may define just as the difference 

/o'"^'(i) = fo^^t) - fl''%). (49) 

We can define the isolated ^"0(1430) resonance contribution to the B Ktt^tt^ decay amplitudes 
by replacing the scalar form factor f^'^{t) by /o°''^(0- For instance, for the B~ K~'K~^Tr~ case, 
this substitution can be done in Eq. (10). 

3. Vector form factor 

We perform a construction using the same method as in the case of the scalar form factor. Here 
the main points involved in this construction are given below, a more detailed discussion can be 
found in Ref. [28]. For Ktt scattering in P-wave at medium energies (m^^ < 2.4 GeV), we assume 
the dominance of two inelastic channels, K*tt and Kp. We treat vector mesons as stable particles 
in our unitary equations. For the three coupled channels Kir (i = 1), K*tt {i = 2) and Kp {i = 3) 
three form factors Hi{t) enter the calculation. The first one was defined in Eq. (7): 

H,{t) ^ V2fr^\t) = fr^\t). (50) 

The other two are defined by the following matrix elements: 

{K*^(j>v,\)\u^^s\'k\p^)) = e^,^f,s*''{X)p^p^,H2{t), (51) 

{p^'ipv, X)\uj^s\K-{pk)) = -e^,uape*''{X)p^p^KHz{t). (52) 

In the above equations is the polarization four-vector of the K* or the p meson, p denote 
four-momenta of mesons and is the completely antisymmetric tensor. The components H2 

and i?3 have dimension mass~^ while Hi is dimensionless. As for the 5— wave case each form 
factor satisfies an unsubtracted dispersion relation. We have now to determine six independent 
matrix elements of the 3x3 P-wave T-matrix, T^, from fits to the experimental data. The most 
complete data exist for the elastic channel K-k —> Kn [23]. Some information is also available on 
the inelastic amplitudes Ktt K*Tr and Ktt Kp in the regions of the resonances if* (1410) 
and if* (1680) [22, 45]. As in the case of the S-wave, the if -matrix method is used to enforce 
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three-channel unitarity. The unitarity equations obeyed by the three form factors Hi{t) can be 
written as 

3 



(53) 



where r, Q and E are diagonal matrices with r = diag[l, \/t, \/t\, E = 2Q/\/t and 

Q = diag[e{t - ti)qK^{t), e{t - t3)qK*At), 0{t - U)qKp{t)] (54) 

with t3 = {niK* + "1^)^ and t/^ = {rriK + rnpY . The normalization of the T^-matrix here is such 
that its relation with the 5'-matrix is 



(55) 



In the asymptotic region, mKir ^ 2.4 GeV, we impose a smooth interpolation for the scattering 
amplitudes with the index = 4 [see Eq. (32)]. In order to enforce exact unitarity in the inter- 
polation region we write the S'-matrix in exponential form, = exp(2i/i), where h is real and 
symmetric matrix and we interpolate the non-diagonal elements of h to zero and the diagonal ones 
to multiples of vr. Illustration of this interpolation of the eigenphases of the S'^-matrix is given in 
the Fig. 4 of Ref. [28]. 
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FIG. 3: Modulus (left panel) and phase (right panel) of the strange vector form factor ^ (t) obtained 
by solving a three-channel Muskhelishvili-Onines equations system. Dependence on the symmetry breaking 
parameter a is illustrated [see Eq. (59)]. 



For = 4, in order to solve the system of integral equations we must impose four conditions 
on the form factors. We use the three values of Hi at t = and, as the fourth constraint, we 
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implement an asymptotic condition for t — oo on the form factor Hi. According to Refs. [20, 55] 
and ignoring flavor symmetry breaking, one should have, 

We do not attempt to reproduce the logarithmic running of ag and actually implement Eq. (56) 
with the constant value ag = 0.2. At t = 0, one has i?i(0) = i^i(O) and we will use the value given 
in Eq. (38). The values of i?2(0) and -^3(0) are not known as precisely. In the SU{3) chiral limit, 
flavor symmetry leads to the following relation between the charged current matrix element and 

the electromagnetic one, Jgj\^ 

{K* + \u^^s\7r') = ^{p+\f^^\^+), (57) 

which allows one to relate i?2(0) and -ff3(0) to the radiative decay width of the p+, yielding 

H2{0) = -i?3(0) = (1.54 ± 0.08) GeV-^ (58) 

The relative sign is determined by vector meson dominance arguments. We have studied the influ- 
ence of flavor symmetry breaking to first order in the quark masses. There are three independent 
symmetry breaking parameters and two of them can be determined from experiment (see [28]). As 
a consequence, one can express H2 and H3 in terms of the third, unknown parameter a, as follows: 

i?2(0) = (1.41 ±0.09 -65.4a) GeV'^, HsiO) = (-1.34 ± 0.07 - 65.4a) GeV-^ (59) 

The magnitude of a is expected to be a few times 10~^. We have estimated a from the sum R of 
the decay rates of the r" into K-jt^ and i?°7r" R{t Kiri/r) = (13.5 ± 0.5) x IQ-^ [52] which 
gives a = (— 7.OI20) ^ 10~^. Results for the vector form factor are displayed in Fig. 3. Its phase 
shows a sharp transition as a function of the parameter a close to a = — 7 x 10~^. It goes from a 
regime where its value is Sir at infinity to one where it is tt displaying a sharp drop. The if* (1680) 
resonance appears to be suppressed but the properties of the form factor in this energy region 
depend significantly on the iS-matrix interpolation parameters in the region rnKv > 2.4 GeV. In 
our application to B decays wc will use the form factor in the region Vi < 1.8 GeV, where the 
sensitivity to the asymptotic interpolation is small. 

4. -K'*(892) part of the vector form factor 

As before, the starting point is the discontinuity equation satisfied by fi'^{t) across the elastic 
unitarity cut, which reads 

/f '^(t + ie) - ff^'it - ie) = -2aKAt + i^hxAt + i^)T[i{t + ie)ff''{t - ie), (60) 
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and from this we deduce the expression of fi^{t) on the second Riemann sheet 



frit) 



(61) 



The position of a resonance in the complex plane, t^"'*^, corresponds to a zero of the denominator 
function in Eq. (61) and this allows one to isolate a pole in ff\t), 



(62) 



f3{t- tf') ■ 

In the case of the i^*(892) we obtain, based on our fit, the following values for the pole parameters 



^poie ^ (0.7982 - i 0.0504) GeV^ -y/tf = (0.8939 - i 0.0282) GeV, 

P = (-1.8874 + i 9.5726) GeV^ /f ^(^r'') = 0-8244 - i 9.0784. (63) 

The i^*(892) resonance being very narrow, we expect the pole component to strongly dominate the 
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FIG. 4: Modulus of the vector form factor {t) compared with the K*{892) complex pole component 



form factor below one GeV. This is illustrated in Fig. 4 which shows that, indeed, the background 
component is very small in that case. As before we can define the B decay amplitude B K*{892)Tr 
from the three-body amplitude B Kirir by substituting fi'^{t) by fi°^'^{t) in the relevant 
formulas, for example in Eqs. (11) and (15). 
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V. FIT, RESULTS AND DISCUSSION 



A. Fitting procedure 

In this paper we use a fitting procedure similar to tliat described in Ref. [19]. We perform a ^ 
fit on experimental data from Belle [5, 6, 8, 9] and BaBar [7, 10-12] collaborations. We use six mKix 
and five cos Qn distributions where the background is subtracted. These data are extracted from 
the figures of the first seven papers just cited. We also exploit the four branching fractions for the 
ir*(892)7r, the three CP asymmetries for if*(892)7r and the three CP asymmetries for K^(1430)7r 
calculated by experimentalists in their data analyzes [6, 9, 11, 12]. Branching fractions are necessary 
to determine the absolute size of decay amplitudes. However, the branching fractions for the B — >■ 
-fCo(1430)7r are not well determined due to the large width of the ^"0(1430) resonance. Therefore 
we use here only the well measured branching fractions for the B —>■ K*{892)it, the i^'*(892) being a 
narrow resonance. The phases of the decay amplitudes can be constrained by the phase difference, 
A$o, between the decay amplitudes of B^ ii'*+(892)7r- and B^ K*-(892)7r+. Here we use 
the preliminary result of Ref. [10] . The total reads 



where the coefficient w is introduced in order to increase the weight of the branching fractions, CP 
asymmetries and the phase difference, which form a significantly smaller data set than the mKn 
and cos 9h distributions. We have verified that varying w between 5 and 20 leads to very similar 
fits. In this analysis, to perform our best fit, we choose w = 10. 
The for ^ given distribution with n bins is defined by 



where Yg^pixi) and AYgxp(xi) are the number of experimental events and associated error in each 
bin Xi. Here x denotes either mK-jrOi cos 9h- Integration of the differential distributions dB{x)/dx 
(see Eqs. (23) and (26)) over the bin width [xj_i, Xj] yields the theoretical number of events Yfh{xi), 



Our theoretical distributions are normalized to the number of experimental events in the analyzed 
range from xq to Xn with 



xL = Xmx. + ^IosOh + ^ (XSR + XAcp + Xa$o) 



(64) 




(65) 




(66) 



n 



N = 




(67) 
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Altogether in our fit we include 308 bins in the uiktv and cos 9h distributions, four branching 
fractions, six CP asymmetries and one phase difference. 

In our minimization, the branching fractions and asymmetries are calculated in limited mKn 
regions. For the K* (892) we choose a mx7r range from 0.82 to 0.97 GeV and for the JCg (1430) one 
from 1.0 to 1.76 GeV. These ranges have been used in Refs. [5, 8, 9] to obtain the helicity angle 
distributions in the regions where the resonances dominate. The experimental branching fractions, 
which we use in our fit, are calculated in the above mi<:7r regions from the models presented in the 
experimental analyzes (see in particular Refs. [5, 9, 11]). In our analysis we have excluded two 
experimental points in the itiktv distributions and three in the cos 9h ones. As will be seen below 
these points lie significantly far from the general trend of the data. 

TABLE II: Phcuomcuological pavamctci-s of the decay amplitudes (see e.g. Eqs. (10). (11) (14) and (15)). 

Real part Imaginary part 

-0.402 ± 0.244 -3.641 ± 0.054 

cl +0.015 ±0.003 +0.033 ±0.004 

-0.051 ±0.153 -0.161 ±0.184 

eg +0.075 ± 0.009 -0.033 ± 0.007 



B. Results 

In Table II we give the values of the phenomenological c^'q parameters and their errors obtained 
from our best fit. These parameters enter our amplitudes defined in Sec. II. The large values of the 
C4 g coefficients should not be directly compared to those of C4 g, since in the amplitudes the latter 
are multiplied by Ac and the former by with |Ac| c± 50|Au|. The results of our fits, with x^= 541 
for 308 experimental Ktt effective mass and helicity angle distribution points, and x^=9.3 for ten 
experimental branching ratios and asymmetries, are presented in Tables III to V and in figures 5 to 
11. For the phase difference A(pQ our fit gives — 199°±6° to be compared with (— 164±24±12±15)° 
found in the experimental analysis of Ref. [10]. In the calculation of distributions wc take into 
account all the tttt effective mass cuts around the D, J /'^ and ^(25") meson masses. As described 
in Refs. [5-11], these cuts are introduced in the experimental analyzes in order to eliminate the 
decay contributions from these resonances. Our results, shown as histograms, take into account 
the above cuts. In all figures presented in this section, the dashed and dotted lines describe the S- 
and P-wave contributions, respectively. In these latter cases the cuts are not taken into account. 



25 



150 h 

> 

^100 



CO 

> 50 

0) 







1 1 1 1 1 

; 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

B* ^ 7T^ 7T* - 


f 


' T 







0.6 0.8 1.0 1.2 1.4 1.6 



FIG. 5: The Jf^^Tr^F effective mass distributions in the — > K^tv^tv^ decays from the fit to the experi- 
mental data as described in Sec. V A. Data points are from Ref. [6] . The dashed fine represents the 5- wave 
contribution of our model, the dotted line that of the P-wave and the histogram corresponds to the coherent 
sum of the S- and P-wave contributions. 



Note that, following the experimental procedure, for the rriKn plots the requirement is made that 
mjr+TT- is greater than 1.5 GeV for the Belle data and 2.0 GeV for the Babar results. Thus, 
the contributions unrelated to Ktt rescattering and arising, for example, from B p{770)K and 
B fo{980)K decays, are removed from the data samples. In our model such contributions are 
omitted. 



1. The Ktt effective mass and helicity angle distributions in decays 

The K^n^ effective mass distribution for the B^ K^tt^tt^ decays, obtained from our fit, is 
compared to the experimental distributions of Belle [6] and BaBar Collaborations [11] in Figs. 5 and 
6, respectively. The mass distributions are averaged over charge conjugate states. In both cases, 
our model describes quite well the m^jr distributions in the (892) and (1430) regions. It also 
depicts quite well the sizable enhancement below 1 GeV related to the Kq{800) state, often called 
K [52, 56]. In our amplitude, its contribution is present in the relatively large background found in 
the modulus of the scalar form factor at low mK-ir (see Eq. (49) and Fig. 2). If one approximates 
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FIG. 6: As in Fig. 5 but for the data points from Ref. [11]. 



the quasi-two body B {Ktt)s vr by the two-body B Kq{1430) tt amphtudes, then is fixed 
at the resonance mass "T'x*(i43o) third hnes of Eqs. (10) and (14). With this replacement, 

one cannot reproduce the low mKn distributions below about 1 GeV, the g^-term contribution 
becoming much too large. This justifies the form of our three-body approach, within the QCD 
factorization framework, to these decays. The origin of this term is given below Eq. (10). 

In the isobar model, used in experimental analyzes, the above dependence is approximated 
by one fitted constant parameter. For the description of the wide ifQ(1430) resonance this may 
not be a very good approximation since the q^ varies by a factor of eight from the Ktt threshold 
to the rriK-K limit of about 1.8 GeV, close to the sum of the mx*(i43o) mass and its width. 

The results of the fit to the cos 9h distributions of the Belle Collaboration [5] around the K*(892) 
and -fCg(1430) resonances are shown in Figs. 7a and 7b. Here we show the contributions of the S 
and P waves as their interference given by Eq. (28) . In the fit corresponding to Fig. 7a, we have 
excluded one bin at cos Oh = —0.95. This bin is not related to any cut and its value is almost 
twice as large as the value of the total for this distribution. In Fig. 7a, the P-wave contribution 
dominates, those of the S wave and of the interference term being rather small. On the contrary, 
and as expected, it can be seen in Fig. 7b that the S-wave contribution is much larger than that 
of the P wave. Contribution of the interference term leads to a visible left-right asymmetry. The 
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FIG. 7: Helicity angle distributions for B"^ — > K'^tt^tt'^ decays calculated from the averaged double dif- 
ferential distribution integrated over mK±nT mass from 0.82 to 0.97 GeV in the K*{892) case a) and from 
1.0 to 1.76 GeV in the /Cq (1430) one b). Data points are from Ref. [5]. Dashed lines represent the S-wa.ve 
contribution of our model, dotted lines that of the P-wave, the dot-dashed that of the interference term. 
The histograms correspond to the sum of these three contributions. 



minima in the histograms of Figs. 7b and 10b at cos6h — —0.15 and +0.2 correspond to the cuts 
related to background events of J/* and ^'(25'), respectively. There is also a cut in Figs. 7a and 7b 
at cos 6h — —0.75 corresponding to the D meson. 



2. The Ktz effective mass and helicity angle distributions for B° —>■ K^tt'^tt and 
B° K°n-TT+ decays 

The results of the fits to the rriKn distributions for the — > K^7r~7r+ and 5° — Kgir'^Tr" 
decays of the Belle Collaboration [9] are shown in Figs. 8a and 8b, respectively. In the fit to the 
B° K'^tt-tt^ distribution, we have excluded the two bins at mKir = 1538 and 1588 MeV. They 
have quite small errors and negative numbers of events. In the mi^: 7^ range around the i^*(892), 
the histogram of the model has less events for case than for B^. One then expects a negative 
CP asymmetry, which is confirmed by the value given in Table IV. It is in agreement with the 
result of the experimental analysis. 

Our rriKTv distributions for the same neutral B decays are compared to those of the BaBar 
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FIG. 8: As in Fig. 5 but for B° K^tt~tt+ decays a), for B° K'^tt+tt~ ones b) and for the data of 
Ref. [9]. 
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FIG. 9: As in Fig. 5 but for 5° K^-k 7r+ decays in a) and for 5° Kg-n^ir ones in b), and for the 
data of Ref. [10]. 
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FIG. 10: As in Fig. 7 but for B° ^ K^Tr^TT+ and B° K^n+n' averaged distributions, and for the Belle 
data of Ref. [9] in a) and that of Ref. [8] in b) . 
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FIG. 11: As in Fig. 7 a) but for B° — * KgW tt"*" and B° KgTT~^Tr averaged distributions, and for the 
BaBar data of Ref. [7]. 
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TABLE III: Branching fractions for the B K*{892)Tr decays averaged over charge conjugate reactions in 
units of 10~^ . In the second column, giving the experimental branching ratios, the 2/3 factor arises from 
isospin symmetry. The values of the model calculated by the integration on niKn from 0.82 to 0.97 GeV 
are compared to the corresponding Belle and BaBar results given in the fourth column. Model errors stem 
from the phenomenological parameter uncertainties obtained through the minimization procedure. The last 
column corresponds to the model without phenomenological parameters. 

Decay mode S^xp j^gf S^xp^Q 0.97) model model [cf = 0] 

B- ^[K*°{892)^ K-TT+jTT- 6.45 ±0.71 [6] 5.35 ± 0.59 5.73 ±0.14 1.42 

7.20 ±0.90 [11] 5.98 ±0.75 

B° ^ [K*-{S92) ^ K°TT-] TT+ 5.60 ±0.93 [9] 4.65 ±0.77 5.42 ±0.16 1.09 

^(11.7 ±1.30) [12] 6.47 ±0.72 

Collaboration [10] in Fig. 9. As previously, the if* (892) and ^^0(1430) are well described by our 
model. Here, the width of the bins is larger than that of Belle in Fig. 8 which explains why 
the maximum of the P-wave contribution is above the experimental points close to the K*{892) 
position. 

The averaged Kg-K^n^ and i^^vr+vr" helicity angle distributions for the mKn 

regions of the K*(892)and i^Tg (1430) are compared in Fig. 10 to the Belle data [8, 9]. In the fit 
to the distribution shown in Fig. 10b we have excluded two bins at cos 6h equal to ±0.95. These 
two data lie rather outside the general trend of the distribution. As in the charged B decays, the 
P wave dominates the K*{892) region and the S wave the Kq(1430) one. In both cases the S-P 
interference is rather small. 

Figure 11 shows our helicity angle distribution fitted to the BaBar experimental data [7] which 
results from the integration of the double differential distribution over uik-k from 0.776 to 1.01 
GcV. The contributions of the S and P wave and of their interference are similar to those observed 
in Figs. 7a and 10a. 

3. Branching fractions and CP asymmetries 

In Table III, our branching fractions for the B —> iC*(892)7r decays are compared to the cor- 
responding experimental Belle and BaBar values. As already mentioned in Sec. VA these are 
obtained from integration on uikit from 0.82 to 0.97 GeV of the double differential branching frac- 
tions of our model and those of the experimental analyzes. The theoretical errors are calculated 
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using covariance matrix elements and the corresponding derivatives of the branching fractions over 
all fitted parameters. We did not include other uncertainties entering our amplitudes, so our theo- 
retical errors are underestimated. Our model branching fractions for B i^*(892)7r decays agree 
quite well with the corresponding experimental ones within their errors. Using the pole part of 
our P-wave form factor (see Eqs. (62) and (63)) and integrating over the full range, we obtain for 
the charged averaged branching fraction 6.95 x 10"^. This value compares quite well with those 
of Belle and BaBar given in Table III. Let us stress however, that when the phenomenological 
parameters (i = 4, 6 and p = u,c), are set to zero, the theoretical branching fractions are too 
small by a factor 4 or 5 (see the last column of Table III). This indicates that QCDF P-wave 
amplitudes are too small by a factor of about 2. 

In Fig. 12a we present the reduced complex P-wave amplitudes, A^p*^, which are given by the 
expression between the curly brackets in Eqs. (11) and (15) for B~ and P° decays, respectively. The 
reduced S'-wave amplitudes, M^g"^, for and B^ decays, shown in Fig. 12b, arc defined similarly 
from Eqs. (10) and (14). The corresponding B^ and B^ amplitudes can then be obtained through 
the conjugation A„ A*, Ac — > A*. The dashed arrows correspond to the reduced amplitudes 
before the fit (c^ = 0) while the solid ones to the result of the fit. In Fig. 12a, the B~ and P+ 
reduced amplitudes without phenomenological parameters are degenerate due to the dominance of 
the almost real Ac term over the A„ one. Moreover, the fact that the B^ and P° reduced amplitudes 
for = have almost opposite imaginary parts, comes from the presence of the tree term A^ ai 
with a large real part of the ai coefficient close to 1 (see Table I). 

TABLE IV: Direct CP asymmetries averaged over charge conjugate reactions. The values of the model, 
calculated over the niKTr range from 0.82 to 0.97 GeV for the Kir P-wave and from 1.0 to 1.76 GeV for 
the 5-wave, are compared to the Belle and BaBar results. Concerning the errors of the model and the last 
column, see the caption in Table III. 

Decay mode exp. (%) Ref. model (%) model (%) [c^ = 0] 

B- ^[K*°{892)^ K-TT+jir- -14.9 ±6.8 [6] -2.5 ±1.3 1.4 

3.2 ±5.4 [11] 
B- ^ [K^{U30) K-TT+j TT- 7.6 ± 4.6 [6] 

B- {K-7:+)s n- 3.2 ± 4.6 [11] 5.4±1.0 0.2 

B^ [K*^{892) ^ K^n-] TT+ -14 ±12 [12] -19.6 ±3.0 6.1 

B° ^ {K°'n:-)s TT+ 17 ±26 [12] -0.2 ±1.3 -1.7 

Our direct CP violating asymmetries are compared with the experimental ones in Table IV. 
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FIG. 12: Complex plane representation of the parts of the amplitudes depending on the CKM matrix 
elements, on the effective QCD coefficients and on the fitted penguin parameters c^. a) P wave, b) S 
wave at mK-R= 1.414 GeV {K^{lA'i{]) mass). Dashed lines: amplitudes with =0, solid lines: results of 
the fit. See text for the exact definition of these reduced amplitudes. 

Their errors are calculated in the same way as for the branching fractions. For B~ — > X*'^(892)7r~ 
decays our asymmetries lie between those of Belle and BaBar. The results for B~ (K~7r~^)s tt~ 
and for {K^7:~)s vr^ decays agree with the experimental values of both collaborations. For 

the /S-wave, the variation with the range of integration is within the experimental errors. For 
instance, if one calculates the asymmetries over the m^Tr range from threshold (0.64 GeV) to 
1.76 GeV, our result for the charged B decays varies from (5.4ibl.0)% to (9.8ibl.l)% and for the 
neutral ones from (-0.2ibl.3)% to (2. Til. 3)% . There is no variation for the P wave. 

C. Discussion 

In the case of the Kq(1430) resonance, it is difficult to extract the quasi- two body branching 
fraction from experimental data due to the presence of a significant background which can mainly be 
attributed to the broad i^Q(800) resonance. Contrary to the P-wavc amplitude, entirely dominated 
by the K*{892) resonance below 1 GeV, the S wave is more complex. This is exemplified in Fig. 2 
by the comparison of the pole part of the scalar form factor to the complete form factor. We remind 
that our S-wave amplitudes are proportional to this form factor as readily seen, for instance, in 
Eqs. (10) and (14). This S-wave complexity results in different parameterizations in Belle and 
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TABLE V: Branching fractions averaged over charge conjugate reactions B — > {KTT)sTr in units of 10~^. The 
second column gives the experimental results. The predictions of our model, calculated by the integration of 
the rriKTr distribution over mxTr from threshold (0.64 GeV) to 1.76 GeV, are compared to the corresponding 
Belle and BaBar results given in the fourth column. In the first two lines, the Belle branching fractions [6, 9], 
calculated with a Breit-Wigner amplitude, are compared to our predictions obtained from the ii'o(1430) 
pole part of the scalar form factor (see Sec. IV B 2). In the last two lines we show the BaBar branching 
fractions [11, 12] for B {K'k)s tt calculated, in their parametrization, with the part of the decay amplitude 
proportional to the K'k S-wa,ve T-matrix. This is compared to the results of our model, where the B — > 
{KTr)sn amplitude corresponds to the part proportional to the scalar form factor (see Sec. IVBl). See 
caption of Table III for the factor of 2/3 in the first column, for the errors of the model and for the last 
column. 

Decay mode ^^^P Ref. ^*^^P(0.64, 1.76) model model [cf = 0] 

B- [K^°{U30) K~TT+] TT^ 32.0 ±3.0 [6] 27.0 ±2.5 11.6 ±0.6 6.1 

i?° -> [i^o~(1430) ^ i?°7r-] 7r+ 30.8 ± 4.0 [9] 26.0 ±3.4 11.1 ±0.5 5.7 

B- ^ {K-7r+)s TT- 24.5 ±5.0 [11] 22.5 ± 4.6 16.5 ± 0.8 7.5 

B° ^ {K'^n-)s TT+ ^(28.2 ±7.5) [12] 17.3 ±4.6 15.8 ±0.7 7.1 



BaBar analyzes. The Belle group uses a Breit-Wigner amplitude to represent the Ktt S-w&ve 
interaction. They have furthermore a large contribution from a nonresonant part. The BaBar 
Collaboration has introduced a term proportional to the Ktt S-wave T-matrix, and used the LASS 
parametrization of the latter. It consists of an effective range nonresonant component plus a 
iiTp (1430) Breit-Wigner term. Since this parametrization is fitted to experimental Kir scattering 
data (in the range 0.8 ^ rriKn ^ 1-53 GeV), this method provides an improved treatment of the 
final state interaction as compared to the Belle parametrization (sec e.g. Ref. [57]). However, 
factorization implies that the B decay amplitude should involve the scalar form factor rather than 
simply the T-matrix. Note that the T-matrix and the associated form factor have the same phase 
in the elastic region (m^Tr ^ 1-45 GeV), but not the same modulus. Above the Ktt elastic region 
both the phase and the modulus are different. 

Our model, based on factorization, allows us to calculate, in an unambiguous way, the branching 
fractions making use of either the complete Kir S'-wave contribution or of the ii'Q(1430) resonance 
only, described here as the pole position of the scalar form factor on the second Riemann sheet (cf. 
Sec. IV B 2). Our predictions, using the rnxn range from threshold (0.64 GeV) to 1.76 GeV, are 
shown in the fifth column of Table V. In the two first lines we use the pole-part contribution of the 
scalar form factor whereas in the last two lines the full scalar form factor contributes. Our values 
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have to be compared with those of the fourth column calculated by us using the experimental 
parameterizations with our range of integration. The values of the experimental analyzes are given 
in the second column. They correspond to integration over the full mKn range. In our model, if 
one integrates also over the full range using the pole part of the scalar form factor (see Eqs. (45), 
(46) and (47)), the averaged charge branching ratio for the decays is 12.7 x 10~^. 

For both charged and neutral B A'g(1430)7r decays, the predictions of our model are smaller 
than the corresponding Belle results [6, 9] by a factor of 2.3. It is worthwhile to mention that 
the Belle Collaboration has found two solutions in their Dalitz-plot analyzes. For example, in 
the Table IV of the second paper of Ref. [6], the value of the solution 2 is smaller by a factor 
of about 5 than that of the retained solution 1. In the case of their solution 1, there is a strong 
negative interference between the resonant Kq(1430) contribution and the nonresonant term. In 
their B^{B^) — > X^tt^tt^ analysis [12] the Babar Collaboration has found four degenerate solutions 
and the quoted errors of the final result (see their Tables IV and V) include the spread of these 
four solutions. 

In Fig. 13, we compare the m^Tr distributions of the averaged B^ {K'^Tr^)s tt^ decays 
corresponding to our model and to the BaBar parametrization calculated with the central values of 
their parameters. As mentioned above, our B — ^ {K'k)s tt amplitude is proportional to the strange 
Kit scalar form factor (see Eqs. (10), (12) and (13)) but in BaBar's case it is the part proportional 
to the K-K 5- wave T-matrix. In Fig. 13a, we show the resonant /^q (1430) contribution (dashed line) 
of our model together with the background part (dotted line) and the interference term (dashed- 
dotted line) between the resonant and background contributions. In Fig. 13b, the corresponding 
three components for the BaBar parametrization [11] are shown, their eff'ective range term (dashed- 
dotted) line playing the role of the background. In our case the interference term is positive, its 
contribution being close to 19%, that of the resonance about 70% and that of the background part 
around 11 %. This can be compared with the BaBar parametrization which, for the range vtik-k 
between 0.64 and 1.76 GeV, gives —25 % for their negative interference contribution, 78 % for 
their resonance part and 47 % for the nonresonant effective range term. These numbers are very 
close to the corresponding values, —26 %, 81 % and 45 % given in Ref. [11] obtained integrating 
over the full range of rriKn- Although these effective mass spectra are significantly different, our 
integrated value for the B~ — > {K~Tr~^)s tt" branching fraction (16.5 ± 0.8) x 10~^ is within one 
standard deviation with respect to the experimental BaBar result (22.5 it 4.6) x 10~^. For the 
neutral B decays the comparison is even better: we obtain (15.8 ± 0.7) x 10"^ and the BaBar 
result, recalculated for the iuktv range from 0.64 to 1.76 GeV, is (17.3 it 4.6) x 10~^ (see last line 
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of Table V). We then suggest below [see Eq. (68)] a parameterization (based on our amplitude) 
proportional to the Kir scalar form factor and which could be used, instead of a parameterization 
proportional to the Kir S-w&ve T-matrix, in experimental analyzes of B — ^ Kwrr decays. 

Our theoretical QCDF predictions (c^ = 0), shown in the last column of Table V, give too 
low branching fractions for all B decays into {Ktt)s vr or Kq{1430) vr by a factor close to 2. 
Figure 12b illustrates the influence of the phenomenological parameters on the theoretical reduced 
S-wave amplitudes. The modulus of the amplitudes increases by a factor of about \/2 and there 
is also an important phase change. The fact that, at the Ktt mass equal to the Kq(1430) mass, 
the magnitudes and phases of all reduced S-wave amplitudes (without and with ) are similar, 
comes from the dominance of the term proportional to Ac and from the smallness of the a^''^{S) 
coefficients (see, for instance, Eqs. (14) and Table I). 
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FIG. 13: Comparison of the different coniponents of the averaged itik-k distributions of the B"^ — > 
{K^iT^)sTr^ decays: a) our model, b) BaBar's LASS parametrization [11]. In this calculation our amplitude 
is proportional to the scalar Kw form factor but that of BaBar is the part proportional to the S-wave Kw 
T-matrix. The dashed lines correspond to the resonant Kq{14:30) contributions, the dotted-dashed lines to 
the background, dotted lines to the interference and the solid lines to their sum. 

As was discussed in Sec. IV B 1 and shown Fig. 1, the S-wave form factor depends on the value 
of the fx/ f-K ratio. We found that the form factors corresponding to fxl fix = 1-193 (form factor of 
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the present work) and 1.183 lead to fits of comparable good quality, with however a slightly better 
for the form factor calculated with /x//7r= 1-183. Use of the form factor with a ratio equal to 
1.203 gives a poorer fit. 

Effective mass and helicity distributions together with branching ratio data allow to determine 
mainly the moduli of the decay amplitudes. Their phases can be constrained by measurements of 
direct CP asymmetries and by time dependent Dalitz-plot analyzes. In our fit, besides asymmetries, 
we use the preliminary value, obtained by the BaBar group [10], of the phase difference A$o 
between the and decay amplitudes to ii'*(892)7r. Had we not imposed this constraint, we 
would obtain four different solutions with equivalent and with almost unchanged moduli for S- 
and P-wave amplitudes but with different phases. In a just published analysis of B^ Kgir'^Tr^ 
decays the Belle Collaboration [58] has proposed values for A$o- We have checked that we can, 
reasonably well, reproduce the value of their solution 2, viz. (14.6;^!^:^ ± 11 ± 17.6)° with a global 
fit of the same quality as the present one. Comparing the results of the two fits, we found indeed 
that the S- and P-wave amplitudes had basically unchanged moduli but modified phases. Note 
the large difference between the preliminary result of Ref. [10] , A$o = (—164 ± 24 it 12 it 15)° and 
the above value of Belle. Our Acp value for P° {K^t^~)p 7r+ given in Table IV, (-19.6 ± 3.0)%, 
agrees well with that of the solution 2 of this Belle analysis, (—20 it 11 it 5 it 5)%. 

As just mentioned above and as found in their previous analysis [6, 9], the Belle Collaboration, 
in this latest time-dependent analysis [58], has retained 2 solutions which is consequence of the 
interplay between the quite broad i^Q(1430) resonance and their phenomenological nonresonant 
background. Their solution 1 has a large i^Q(1430)7r fit fraction with a sizable negative resonant- 
nonresonant interference term while their solution 2 is characterized by a fit fraction smaller by 
a factor of 3.5 and a small interference term, as we found in our model. Using the i^T^Tr+Tr" 
charmless total branching fraction (47.5 ± 2.4 it 3.7) x 10~^, as given in Table III of their previous 
publication [9], one obtains a branching fraction of 8.3 xlO~^ for their solution 2, value close to 
our result cited in Table V. 

VI. SUMMARY AND CONCLUDING REMARKS 

With this analysis of S- and P-wave pion-kaon interactions in S ^ K-ktt decays, we have 
extended and completed previous studies on final-state interactions in these three-body de- 
cays [17, 19]. We have concentrated on the scalar i^Tg (1430) and vector K*{SSi2) resonances following 
the logic of these previous works that treated the tttt interactions in S- and P-waves as well as their 
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interferences. In Sec. II the weak decay amplitudes were again derived in the QCD factorization 
approach [13, 32], which express them as a product of two currents multiphed by a sum of effective 
coefficients which includes non-factorizable contributions. These coefHcients, representing pertur- 
bative QCD leading order amplitudes and next-to-leading order vertex and penguin corrections, 
were studied and given in Sec. III. The contribution of these higher order terms is not sufficient 
to obtain a good description of data. Therefore, we have introduced phenomenological parameters 
which can simulate on one hand long-distance charming penguin amplitudes [30] and on the other 
hand hard-spectator scattering and weak annihilation contributions [13]. These phenomenological 
amplitudes could receive also, through 6 to s quark transitions involved here, some new-physics 
contributions. 

The different models for the matrix element of the first current, expressing the B to tt transition 
in terms of the scalar and vector transition form factors have been briefly reviewed in Sec. IV A. 
The creation of a pion-kaon pair in an S- or P-wave from vacuum is mediated by the second 
current, and accordingly described by a Kit scalar and a vector form factor. These control the 
dependence of the decay amplitude as a function of the Kir invariant mass, because the Bn form 
factors are nearly constant in the region considered. In Sec. IV B 1 the scalar Ktt form factor was 
calculated along similar lines as in Ref. [26] and the extension to the case of a vector form factor 
was developed in Sec. IV B 3. 

We treat both the S- and P-wave on the same footing, namely relating the form factors using 
their analyticity and unitarity relations to pion-kaon scattering properties known from experiments. 
A simplified, but realistic treatment of inelasticity is also implemented. In the determination of 
these form factors we also use chiral symmetry and QCD constraints at low and high energies 
respectively. As a byproduct of the scalar form factor study we predict for the modulus of the 
1^0(1430) decay constant a value of 32 ± 5 MeV. Our theoretical amplitudes go beyond the usual 
two-body approach applied to decays such as B ^ K*Tr and correctly accounts for the Ktt final- 
state interaction both on and away from resonance peaks. A nonresonant background can be 
isolated from the resonant one in our model as illustrated in Sec. IV B 2, and compared to those 
introduced by Belle as well as BaBar collaborations in their amplitude parametrization [6, 9, 11, 12]. 
A comparison between our resonant, nonresonant and interference term splitting and that of BaBar 
was presented in Sec. VC. 

Furthermore, our model correctly reproduces the enhancement of the decay distributions in the 
low-mass region as observed in Figs. 5, 6, 8 and 9. This enhancement may be attributed to the 
broad Kq{800) resonance which is present in the T-matrix that we use. As can be seen from Figs. 2 
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and 13 the (800) is responsible, in our model, for the behavior of our S'-wave amplitude for rriKTr 
from threshold to about 1.2 GeV. 

Our theoretical QCDF amplitude predicts branching fractions for the B K*{892) tt and 
B —>■ Kq(1439) tt decays too small by factors of about 5 and 2, respectively. The inclusion of 
four complex phenomenological penguin parameters allows us to have a realistic model. These 
parameters, common for B^ , B^ , B^ and B^ decays, have been fitted to numerous experimental 
data, which includes 308 data for the Kn effective mass and helicity angle distributions, four 
B K*{892) TT branching fractions, six direct CP violating asymmetries and the phase difference 
between the B^ and B^ decay amplitudes to K*(892) tt. Our model reproduces rather well these 
319 data with total of 551.5 corresponding to a per degree of freedom equal to 1.77. This 
good reproduction of the data makes it possible to predict the B {Kit)s tt and B Kq(1430) tt 
branching fractions. We can obtain, without ambiguities, the pole contribution of the iiTp (1430) 
resonance. This contribution is smaller than the experimental determination by Belle (see Table V) 
and BaBar [12] by factors of 2.3 and 1.4, respectively. The determination of these branching 
fractions, within the isobar model is problematic since the resonance i^Q(1430) is wide and the 
nonresonant part difficult to assess. The non-uniqueness of the parametrization of the nonresonant 
contribution leads to a large systematic uncertainty of the B — ifQ(1430)7r branching fraction as 
seen in the Particle Data Tables [52]. In our approach, with a scalar form factor well constrained 
by theory and experiments other than 5-decays studies, we can describe the data over the rriK-K 
range from threshold to 1.8 GeV. 

Therefore, to diminish ambiguities in data analyzes, we propose to use the following S-wave 
amplitude parametrization for B — > {K'k)s'^ decays: 

Ms{mK.) = ff^ml,) r^ + ciV (68) 

which follows from Eq. (10). Here, cq and c\ are constant complex parameters to be determined 
through the Dalitz-plot analysis of a given 5-meson decay. Upon request, we can provide a numer- 
ical table for the scalar form factor f^'^ {m\,). To calculate the Kq{143>Q) resonance contribution, 
one can replace, once the cq and ci parameters have been determined, f^^ {j^I-k) by its pole part 
/r'*" {m]i,) given in Eqs. (45), (46) and (47). 

The direct CP violating asymmetries and the time dependent CP asymmetries are related to 
the not very well determined angle 7 of the unitary triangle. Our amplitudes are sensitive to 7 via 
their dependence on (see Eqs. (10), (11), (14) and (15)). Precise measurements of the Dalitz 
plot distributions could allow to constrain 7 using our model. 
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APPENDIX A: VERTEX AND PENGUIN CORRECTIONS FOR PS AND PV FINAL 
STATES 

Here we show how we calculate the next-to-leading order vertex and penguin corrections entering 
the effective QCD amplitudes a^ifi) of Eq. (31). We compute their values for PS and PV final 
states, namely KQ(1430)7r and if*(892)7r. Let us first quote some general results from Ref. [13] for 
the vertex correction terms Vi{M): 



Hi = 6, 8, where $m and are the leading-twist and twist-3 distribution amplitudes, respectively, 
of the emitted meson. The integration is over the longitudinal meson-momentum fraction x. As in 
Sec. HI, M stands for the emitted meson that does not include the spectator quark. The functions 
g{x) and h{x) are given in Eq. (38) of Ref. [13]. 

Light cone distribution amplitudes (LCDA) for scalar mesons were derived making use of QCD 
sum rules [33] . We use these distributions to calculate vertex corrections for the case where M is 
a 7^0(1430). The leading twist Gegenbauer expansion for scalar mesons is given by [33] 




(Al) 



for i = 1,4, 10 and 




(A2) 



00 



$5(x) = 6x{l -x) al + J2 «n(/^)Cn/'(2x - 1) 



(A3) 



n=l 
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while the scalar twist-3 amplitude is ^s{x) = 1. Here, the a^{n) are related to the scalar Gegen- 
bauer moments BnifJ-) by virtue of 

with the quark masses rrii ^ m2, as in the case of the i^g(1430), and where nis is the scalar 
meson mass. The normalization condition $5(.T)dT = 1 yields Bq = ji^^ and thus = 1 if the 
small even Gegenbauer moments are neglected. With this, the vertex corrections to order af for 
a ii'^(1430)7r final state with M = if^ are given, for i = 1,4, 10, by the Eq. (4.3) of Ref. [33] and 
byyi(if^) = -6forz = 6,8. 

The leading twist LCDA for a vector meson M = V reads 



^v{x) = 6x(l — x) 



(A5) 



l + ^ar(/x)C3/2(2x-l) 

n=l 

which is given by the same expansion in Gegebauer polynomials Cn^^{2x — 1) as the one for pseu- 
doscalar mesons ^p{x) but with different moments ajf(/x). Thus, the vertex corrections Eq. (Al) 
for a i^*(892)7r final state where M = K*, are given by the Eqs. (47) and (48) of Ref. [32] for 
i = 1,4, 10. Taking into account the twist-3 LCDA, Pn{x) being the usual Legendre polynomials, 

oo 

"n±(/^)-fn+l(2x - 1), (A6) 

n=0 

we obtain the i = 6,8 corrections as 

Vi{K*) = (9 - 6m) af; + (^^ - ztt) a^l , (AT) 

where we have made use of the property ^v{x)dx = 0. 

At order a^, corrections from penguin contractions with the various operators Oi{fx) exist for 
i = 4,6 for QCD penguins and i = 8, 10 for electroweak penguins but not for z = 1. The expressions 
for these contributions can be found in integral form in Ref. [13] for the B PV decay in Eqs. (39) 
to (46). We apply them using the latest results on Gegenbauer moments for the iir*(892) [59]. The 
PS penguin corrections have the same expressions as those for PP final states [32] but one must 
employ the LCDA introduced in Eq. (A3). Nonetheless, since even Gegenbauer moments are 
suppressed, we take into account corrections up to (/x) as for the vertex corrections. 

Finally, the input parameters entering our computation of the a"''^(/i) include the u-, s-, c- and 
6-quark masses, the strong coupling constant ag and the Gegenbauer moments of the leading twist 
and twist-three light cone distribution amplitudes for the scalar ifQ(1430) and vector K*{892) 
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TABLE VI: Next-to-leading 


order vertex and penguin corrections entering a^{mb) [see 


Eq. (31)] for B 


i4:o*(1430) 


77 and B — > if*(892)7r decays. Note that there 


are no penguin corrections for i 


= 1. 




B 


■ (1430)77 


B^K*{892)w 




Vertex 


Penguin 


Vertex 


Penguin 


fll 


0.011 + i 0.063 


0.0 


0.028 + i 0.014 


0.0 




-0.001 - i 0.005 


-0.029 - i 0.019 


-0.002 0.001 


0.004 - i 0.014 


4 




— n (T^T -1- 9 n nfii 

U.UOl t U.UU-L 


-0.002 0.001 


-0.002 -i 0.004 




-0.0004 + i 


-0.003- i 0.014 


0.001 -i 0.001 


-0.007 -i 0.001 


ag 


-0.0004 + i 


-0.006- i 0.004 


0.001 - i 0.001 


0.001 + i 0.011 


"8 


0.0 + i 0.0 


0.0 -i 0.0001 


-0.00001 + i 0.00001 


-0.0 + i 0.0 


as 


0.0 + i 0.0 


0.0- i 0.0 


-0.00001 + i 0.00001 


-0.0 + i 0.0001 


u 
"lO 


0.0006 + i 0.0032 


-0.0006- i 0.0001 


0.0014 + i 0.0007 


0.0002 -i 0.0001 


"10 


0.0006 + i 0.0032 


-0.0007 + i 0.0003 


0.0014 + i 0.0007 


0.0002 - i 0.0 



mesons. We use the scale fi = mjj with as{mh) = 0.223. The corresponding values of the quark 
masses have been given below Eq. (13). In order to calculate the Gegenbauer moments associated 
with the LCDA of the scalar meson, we start from the values at /x = mb/2 from Table X of Ref. [33]. 
For the vector meson we use the recently determined moments af* and af^ at /x = 1 GeV [59]. 
After evolution to the scale /j, = mb, one obtains the following values: ai(Ko) = 5.26, as{KQ) = 
-8.24 for the scalar meson i^o*(1430) and ai{K*) = 0.018, a2{K*) = 0.065, ai±(K*) = 0.026, 
0(2±{K*) = 0.065 for the vector meson i^*(892). In Tabic VI, wc give our results for the ncxt-to- 
leading order vertex and penguin corrections of Eq. (31) and from which are calculated the a^{mb) 
listed in Table I. 

APPENDIX B: DETERMINATION OF THE S-WAVE T-MATRIX ELEMENTS 

Below, we describe the determination of the S-wave T-matrix elements Tn, Tu and T22, the 
channel Ktt being labeled as 1 and Krj' as 2. 

1. Fit above the inelastic threshold 

Precise experimental data on Kir scattering is available [23, 44] and cover approximately the 
range 0.9 ^ ttikw ^ 2.5 GeV. Experiment shows that inelasticity effectively sets in at the rj'K 
threshold and we make the assumption that it is saturated by the ij'K channel. Above the inelastic 
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threshold, the three components of the T-matrix, Tn, Tu and T22 are needed in the unitarity 
equations. A two-channel S-matrix which is unitary and satisfies time reversal invar iance can 
be parametrized, in terms of three observable quantities: two phase-shifts {Skt^ ^Kt)') and one 
inelasticity parameter (rjK-K), 



The relation between the partial-wave T and S matrices is 

Smn + 4i fSnii^lSlM] ' r^„, m, n = 1, 2 with s = {pk + P^f = ml^, (B2) 



-'mn — "mn 



S 

and where qi{s) is the center-of-mass momentum for channel i. A simple way to enforce unitarity is 
to use a iC-matrix type representation of the T-matrix. We take here the following representation 

= K-^ - diag{Ji{s), J2(s)), (B3) 

where K must be real and symmetric and the functions Ji{s), i = 1,2 read: 



, ^ f°° ds' 



with si = {rriK + JTItt)^ and S2 = {rriK + m^/)^. Following the approach of Ref. [26], we use the 
following parametrization for the iT-matrix: 

which includes two resonances and a background term. We determine the parameters by performing 
a fit over an energy range, 1.25 < niKn < 2.5 GeV. We used the experimental data of LASS [23], 

who measured the phase and the modulus of the charged amplitude K^tt~ K^tt^ combined 
with the earlier measurements of Estabrooks et al. [44] of the isospin 3/2 component. The central 
values of the parameters determined from the fit are as follows: 

Ml = 1.454, gi = 0.505, 52 = 1-651, 
M2 = 1.988, hi = 0.784, /12 = 1.144, 

(B6) 

an = 2.371, = 10.060, 022 = -43.946, 

611 = -1.345, 612 = -2.051, 622 = 14.538. 

All the above parameters are in units of GeV except for which are dimensionless and 6^ which 
are in units of GeV~^. The energy cutoff parameter c is not fitted, it is set to c = 1 GeV^. The 
number of data points is 70, the total is 205 and the number of parameters in this fit is 12. 
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Figure 14 shows the comparison of the fit result with the experimental data of Ref . [23] . Using these 
results, when t — > oo, the phases and inelasticity parameters of the S'-matrix satisfy: dKnioo) = 2tt, 
^Kri'ioo) = and rfK-wioo) = 1. 




1.4 1.6 1.8 2 2.2 2.4 

m/tv (GeV) 




1.4 1.6 1.8 2 2.2 2.4 

m/t-TT (GeV) 



FIG. 14: Results of the fC-matrix fit for the modulus, |ao|, and the phase, (/)o, of the K^tt K^tt 
amplitude. 



2. Below the inelastic threshold 

a. Chiral symmetry constraints on T12 

Chiral symmetry constrains scattering amplitudes which involve the r]' meson at low energy if 
one combines chiral symmetry with the large Nc expansion [60]. We will use here a systematic 
expansion scheme based on counting 1/Nc on the same footing as a chiral factor [47, 48] 

(B7) 

At order 5 of this expansion the Lagrangian contains three independent terms 

A = ^ {tr {D.UD^^U^) + tr {x^U + Uh)} - (t>l (B8) 

where [/ is a unitary matrix which contains a nonet of pseudoscalar meson fields (po^-j't's- At this 
order, r] — r]' mixing involves one angle 9 and its value may be determined such that the physical 
r] mass is reproduced. This gives 6 = —5.6° and the prediction for the mass of the t]' meson is too 
large [61] (M^/ ~ 1.6 GeV). This problem is cured by going to the next order of this expansion. 
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The Lagrangian at order 5"^ contains eight independent terms 

Cp = L2tr {D^U^D^UD''U^D''U) + (L3 + 2L2)tr {D^U^ Df'UD^U^ D^U) 

+L5tr {D^U^D^Uix^U + U^x)) + Lstv {x^Ux^U + xU^xU^) 

-iLgtT {F^^D''UD''U^ + F^^D'^U'f D^U) + Liotr {U^ F^^U F'"' ^) 

+hD^^oD^^o + ik2^(l)o tr (x^U - U^x)- (B9) 

This Lagrangian involves the subset of the Gasser-Leutwyler [49] coupUng constants Li which are of 
leading order in Nc (the scale dependence shows up here at order 5^) plus two additional couplings, 
ki and ^2- At order 5^, r] — r]' mixing involves two angles ^0 and Og, 

4>8 = T-(</'r; COS 6*8 + 4>^' sin6'8), 
^8 

^0 = T-{-4>v sin ^'o + (t>n' COS 9o). (BIO) 

The factors Aq and Ag can be expressed in terms of L5, and the angles 9o, 6s can be expressed 
in terms of the chiral couplings L5, Lg and the physical meson masses m^, m^i [47]. Using the 
Lagrangian's (B8) and (B9) a small calculation yields the scattering amplitude Ktt — Krj'. It can 
be written in the following form: 

TKn,Kr]'{s,t, u) = siuOg Ts{s, t,u) + cos9oTo{s, t,u) + sin9Ts{s, t,u) + cos ^Tc(s,i, u), (Bll) 

with 

a/3 

T8{s,t,u) = -—^{-9t + 8rnj^ + ml + Sml,), 
36/^ 

/fi 

Ms,t,u) = ^{2mj, + ml), (B12) 



and 



Ts{s,t,u)=L3-j^[-2{t-2ml){t-m^^,-ml) 

+{s ~ ml - mj^){s — rn^i — m\) + (tx — — m\){u — — m\) ] 



+L8 -^{^ {m\ - ml) {2m\ + ml), (B13) 
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T,{s,t,u) = (L3 + 2L2)-^[(t-2m|)(t-mJ, -m^) 

+{s - ml - mj^){s - m^i - m\) + (u - ml - m\){u - m^, - mff)] 

+-^5 [ ^^{'^'k - '^V) + ^^rn\m1 + 3m|^m^, + 4!m%^ + 6ml,ml - 8m^] 

^ '/tt 

with k2 = k2 — ki/2. Projecting eq. (Bll) on its / = partial wave gives ri2(s). Its value in 
numerical form at the Kit threshold si is at order (S + 6^): 

Tuisi) = 0.32 sine's + 0.28 cos 6*0 + ^2 (-470.5 cos e") 

+L3 (-156.8 cos e + 14.5 sin 6) + L5 (-229.1 cos 6* - 21.5 sin 9) 

+L8 (161.4COS0 - 169.3 sine) + h (0.56 cos 0). (B15) 

The values of the couplings L5, Lg can be determined from the ratio of the decay constants fK/fn 
and the ratio of the quark masses 2ms/{mu + m^) using the 5 expansion up to order (5^. Using, 
for instance, the central values obtained from lattice QCD by the MILC collaboration [62] yields: 
L5 ~ 1.97 X 10-3 and Lg ~ 0.87 x 10"^, while for the mixing angles one obtains 9q c± —18.9° and 
08 — —3.03°. Fitting the 77 and 77' masses in the 5 expansion gives ^2 — 0.12 [47]. Finally, we need 
the values of L2 and L3. In the ordinary chiral expansion, L2{n), L'^i/J-) can be obtained either 
from sum rules based on tttt scattering [49] or based on Kn scattering [63, 64] or from data on 
decay form factors [65, 66]. For illustration, let us adopt the values from [64] and identify L2, 
with L2{fi), Ll^il-i) at = m^. This gives: L2 ~ 1.3 x 10^^ and L3 ~ —4.4 x 10^^ . We can now 
deduce the value of the transition matrix clement T12. At leading order one finds: ri2(si) ~ 0.25, 
while including next-to-leading order corrections one obtains: Ti2(si) ~ 0.15. Clearly, convergence 
is not very fast but we can expect the order of magnitude to be reasonable. This result will serve 
us in the construction of 712 (s) in the unphysical region s < {mx + m^/)^. 

Finally, it is instructive to calculate the predictions for the values of the scalar form factor 
components -Fi(O), ^2(6) in this approach. A small calculation using the Lagrangian gives 



^K-^l' f . 8V2( 



Fi(0) = 1, ^2(0) = 1 { sin^g - -v^v"^^ rn,, ^^ ^^^^^ ^ 

mj^-mi y Sf4 J 

The deviation of -Fi(O) from 1 is proportional to (m^ — m)^ according to the AdemoUo-Gatto 
theorem [67]. In the 6 expansion approach, the deviation shows up at order because it is 
subleading in Nc- The value which we obtain for -P2(0) at order (5^ turns out to be very similar to 
the one obtained in Ref. [26] in a somewhat different approach. The corrections to ^2(0) of order 
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S^, however, have no reason to be particularly small. In fact, the value of ^2(0) which we obtain 
from the solutions of the Muskhelishvili-Omnes equations solutions, using chiral constraints on the 
component Fi{t), is somewhat smaller than in Eq. (B16): -F2(0) ^ 0.52. An analogous result was 
obtained in Ref. [26]. 

b. Determination of Tn and T12 

Below the inelastic threshold, the two components Tn and T12 arc needed in the unitarity 
equations for the scalar form factors. At first, let us look at Tn. At low energies, niK-TT ^ 0.9 
GeV, K-K —>■ Kit scattering (analogously to tttt — tttt scattering) is constrained by Roy-Steiner 
equations [68] which result from combining dispersion relations and crossing symmetry with elastic 
unitarity. We will use the results obtained from a recent re-analysis of such equations [64]. In the 
energy range 0.9 < uikit < 1-25 GeV we also use the fit performed in that reference of the Kir 
elastic phase shift. 

There remains to discuss the Kit Kjf transition matrix clement T12 in the unphysical region. 
For this purpose, we may use a simple method which exploits the fact that the phase, 612, of T12 is 
known for all values of the energy. Indeed, it is identical with the elastic phase below the inelastic 
threshold by Watson's theorem and equal to Sktt + Sxif above because of two-channel unitarity. 
One can then compute the Omnes function: 



and study t/^is) = i^^2 {s)Ti2{s). The function tp{s) has no right-hand cut since Imtp^s) vanishes 
in the range {uik + J^tt)^ < s < 00. Therefore, over a finite interval, we can approximate ipis) by 
a polynomial. In practice, we use a polynomial of degree two and deduce the three parameters of 
the polynomial from the known values of ip{s) at three points, s = si, S2 and S3 = (1.8)^ GeV^. 
The value at the Kir threshold si is (approximately) known from the discussion above using the 
5 expansion, while the values at S2 and S3 are known from the if-matrix fit. The result obtained 
in this manner for the modulus of T12 is displayed in Fig. 15. This completes the determination of 
the three matrix elements Tij{s) in the energy region where they are needed in Eqs. (35) and (36). 




(B17) 



47 




^ (GeV) 

FIG. 15: Absolute value of the T-matrix element Ti2(s). In the region s > {rriK + m.rj') it is obtained from 
the /^-matrix fit, and in the region s < {jtik + it is computed by a polynomial approximation using 
the Omnes function as discussed in the text. 

APPENDIX C: DETERMINATION OF THE DECAY CONSTANTS Jk^ AND fx- IN 
THE COMPLEX POLE APPROACH 

In this appendix, we quote the results for the decay constants fx* and fx* which are asso- 
ciated with matrix elements involving the scalar meson ii'Q(1430) and the vector meson if* (892) 
respectively, in the complex pole approach. For the decay constant fx*, let us use the definition 
proposed by Maltman [69]: 

{0\r^{x)\K*{p)) = fK'^ ml, exp(-zpx) (CI) 

with J^^{x) = dfj_s{x)'y^u{x). We introduce the two-point correlation function associated with this 
current: 

n"^(t) = i J d^xexp(ipx)(0|r[J""(x)(J^")t(0)]|0). (C2) 

Using (CI), the contribution of the ifg(1430) to this correlation function, if it were a stable state, 
would be: 

^K* fx* 

n"*(*)lx*(i43o) = ^2° (^^) 

"■0 

In reahty, the -ftrg(1430) is a resonance and it shows up as a pole of n''*(t) on the second Riemann 
sheet. By analogy with Eq. (C3), we can identify the decay constant fx* from the residue of the 
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pole (which imphes that it is a complex number). As before, in order to define the extension to 
the second sheet we consider the discontinuity of the function 11"* (t) along the real axis: 

n*«(i + ie) - n*«(t - ie) = -^{mj, - mlfaKnis + ie)fo{s + ie)fo{s - ie) (C4) 

for t real and lying in the range between the Kir and the Kr]' thresholds. The factor 3 comes from 
summing over the two possible charge states of the Ktt system. From Eq. (C4), one deduces that 
the extension to the second sheet must be defined as follows: 

167r(l-2(7if^(t)Tfi(0) ^ ^ 

We recognize again here the denominator function D{t) which has a zero at t = to- We can identify 
the residue of the pole at t = to with Eq. (C3), replacing m|^* by (to)^- The following expression 
for Jk* results: 



Numerically, using the preceding results Eqs. (46) and (47) we obtain for the SkI S-k = 1-193, 

fK* = (31.3 + i 7.6) MeV. (C7) 

The result is quasi real and comparable with the value obtained by Maltman [69], /k* = 42.2 MeV. 
Varying fx/ fw one obtains: 

fKS = (36.7 + i 7.5) MeV if ^ = 1.203 and 

Jks = (25.8 + i 9.9) MeV if ^ = 1.183. (C8) 

Let us now present the analogous results for the vector form factor ff'"{t) and the i^*(892) 
resonance. We want now to identify the decay constant fx* associated with the vector meson 
K*{892), which may be defined as 

(Ob7(x)|i^*+(p)) = fK* exp(-ipx) (C9) 

with jf^{x) = s{x)^^u{x). For this purpose, we investigate the correlation function 

= i J d^xe*'^"(0|T(j7(a;)jr (0))|0) = (g^g. - q'g^.u)Mq^) + q,.quM<l')- (CIO) 

The discontinuity along the Kir elastic cut of Ili{q^) reads: 

Ui{t + ie) -Ui(t- ie) = -^^^KAt + ie)qlAt + ie)f^''{t + ie)f^''{t - ie), (Cll) 
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which allows, as before, to obtain the definition on the second sheet 

One then identifies the pole in this expression with the one generated by a stable K* using (C9), 
which, finally, gives fx* in terms of the vector form factor ff^'" 

/ \ 2 

pole\ I i,pole\ fKiT {mPoIc' 



fh = ^ 1 —■ (C13) 

Numerically, inserting the values for t^^^ , fi'^(t^^'^) and (5 (see Eqs. (63)) gives 

fK' ~ (213.9 - i 13.6) MeV. (C14) 

The modulus of this decay decay constant, 214.3, is close to the value fv = 218 MeV we use in 
the P-wave amplitudes (see Eqs. (10) and (15)). 

APPENDIX D: EFFECTIVE DECAY CONSTANTS AND TWO-BODY AMPLITUDES 
B- K*°(1430)n- AND B' K*'^{892)n- 

We shall discuss below the case of K^Q''(1430)7r^ decays. It can be easily generalized to 

decays of other B mesons. For the two body B~ — > ^o°(1430)7r~ decay mode the branching B2S 
can be written in terms of the two body amplitude M2S' 

where |p^*| is the modulus of the ^q(1430) momentum in the B~ rest frame: 



(rriK* + m^) M| - (m^* - m^) 



(D2) 



^0 ' 2Mb- 

ruf^* being the Kq(1430) mass. In the vicinity of mx--K+ = ^R* the two body amplitude M.2S is 
related to the three-body one Ai'^ [Eq. (10)], 

•^5 = ^2S ^%o^K-n+i^K-n+)- (D3) 

The vertex function can be expressed in terms of the scalar form factor /(f- '^'^ (m|-_^+) defined 
in Eq. (7): 

r|-^^-.4-K-..) = j-^ ^4^/0^"'^' (D4) 

° /ifo*(1430) ^K-n+ 
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In this equation /x*(i43o) represents the not very well known ^"0(1430) decay constant (see Ap- 
pendix C and also Ref. [33]). If the mx-Tr+ effective mass is far from the resonance mass m^*o then, 
as one can see from Eq. (10), the relations (D3) and (D4) cannot be used anymore, in particular, 

close to the Ktt threshold and for itik-k ^ "^i?*(i43o)- 

— 

Integration of [Eq. (24)1 over the mx-n+ range from mmin to mmaxj where the 

Kq{143>Q) dominates, gives for the process B~ —>■ {K~Tr~^)sTT~ the branching fraction 



\M 



251 



|2 



/■mmax 

/ dmK-K+ mK-T,+ \'£>i^+\ IPtt-I |r|*o^jf-^+(mK-,r+)l^- (D5) 



4(27r)3M3r^ 

The branching ratios B2S and 833 are simply related by 

Bss = bB2S (D6) 

2 

where b = - 0.93 is the secondary branching fraction for the decay Kq(1430) K^tt'^ [52]. Using 
o 

Eqs. (Dl), (D4) and (D5), one obtains for the modulus square of the effective decay constant 

The knowledge of the scalar form factor allows one to calculate this effective decay constant. 
Integration over the range mmin ^ 'mKn ^ "mmax with mmin = 1 GeV and mmax = 1.76 GeV gives 
1/^/ J = 31 MeV. This value is close to that of the decay constant calculated using the pole 
part of the scalar form factor, 1/)^* (1430) I = 32 MeV [see Eq. (C7)]. This agreement is expected as, 
in the mK-K range (1, 1.76) GeV, the ^g(1430) pole part dominates (sec e.g. Fig. 13). 

The two-body B' i^*°(892)7r- decay amplitude, M^p, can be expressed in terms of the 
three-body one Mp [Eq. (11)], 

Mp = 2 M^p r^.^K-n+i^K-n+)- (D8) 
The vertex function is related to the vector form factor ff '^^(m|._^+) defined in Eq. (7): 

mK*jK* 

where rriK* and fx* are if* (892) mass and the decay constant, respectively. The two-body branch- 
ing fraction for the B" — K*°(892)7r~ decay is 

B2P = \M^p\' „ 'Pf^l! , (DIO) 
51 



where \pk*\ is the modulus of the if* (892) momentum in the B~ rest frame. It can be cal- 
culated from Eq. (D2) replacing m^* by uik*- The three-body branching fraction B^p for 
B~ — K*'^(892)7r~, K*^ — K'''k'^ is obtained by integration of the P-wave part of the effec- 
tive mass distribution [see Eq. (24)] from m^jj^ to m^g^^, covering the range where the K*{892) 
resonance dominates, 



03P = / -; — dm 



K-TT-^ 



\M 



|2 



/ dmx-jj-+ mx-jj-+\p^+\ IPtt-I \^ K*^K--!r+ 



3(27r)3M|r^ 

Now Bsp = {2/3)B2P, where the factor 2/3 is the secondary branching fraction for the decay 
^*''(982) — ifTT"*". As previously for the case of the Kq{1430), taking into account the limited 
range of mx-K between m^j^ and m^^x ^^'^ using Eqs. (D9), (DIO) and (Dll), one obtains the 
modulus square of the effective decay constant f^/ 

-1 /""^Sax 

1/1:^-^1' = V^lJTi ^ / dmK-^+ m^-^+lp^+f |p,-|3|/f""^ (D12) 

^ ' ' min 

Integration from 0.82 to 0.97 GeV, range where the i^*(892) dominates, gives = 194 MeV. 

This value compares well with the decay constant calculated from the pole part of the vector form 
factor, \ fK'* \ = 214 MeV [see Eq. (C14]. A larger range of integration will improve the agreement. 
In the limit of infinite Mb mass and of zero width if* (892), the effective decay constant equals 
fx*- 
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